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ABSTRACT 

A study  is  presented  of  a vortex  formed  near  the  inner  wall  of  a 
rotor  operating  within  a boundary  layer.  Experimental  critical  cavita- 
tion data  and  Laser  Doppler  Anemometer  measurements  obtained  down  tream 
of  the  rotor  plane  are  given  for  many  flow  configurations.  In  addition, 
rotor  inlet  and  outlet  velocity  profiles  were  measured.  Experimental 
results  show  an  influence  of  inlet  velocity  gradient  near  the  inner  wall 
on  the  critical  cavitation  number  of  the  vortex.  Secondary  flow 
vorticity  equations  are  developed  and  employed  to  calculate  the  vorticity 
in  the  blade  passage  near  the  inner  wall  of  a rotor.  In  order  to 
calculate  the  streamwise  vorticity,  an  approximate  method  was  developed 
to  calculate  the  streamlines  through  a rotor.  Cavitation  data  are 
correlated  with  the  calculated  exit  passage  streamwise  vorticity  using  a 
simple  vortex  model.  Results  show  the  importance  of  the  additional 
secondary  vorticity  on  the  exit  streamwise  vorticity  near  the  inner  wall 
of  a rotor.  The  secondary  vorticity  indtices  an  additional  component  of 
tangential  velocity  which  was  found  to  increase  the  flow  turning  near  the 
inner  wall.  This  increased  tangential  velocity  profile  compared  favor- 
ably with  the  measured  vortex  profile  outside  of  the  vortex  core  region. 
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1.  INTRODUCTION 


1 • 1 Ortwin  of  Investigation 

Secondary  flows  generate  additional  streamwise  vorticity  near  the 
Inner  wall  when  a boundary  layer  flow  Is  turned  by  a rotor.  The  apparent 
effect  of  this  additional  vorticity  is  evidenced  by  the  high  cavitation 
numbers  of  the  vortex  formed  along  the  inner  wall  of  a rotor.  One  ex am pi 
of  the  cavitation  associated  with  a vortex  can  be  found  in  the  draft  tube 
of  a Francis  turbine  operating  in  the  part  load  range.  The  cavitation 
number  is  proportional  to  the  square  of  the  streamwise  vorticity  associa- 
ted with  the  vortex.  In  most  cases,  critical  cavitation  numbers  typical 
of  this  vortex  are  often  higher  than  those  associated  with  any  other  type 
of  rotor  cavitation. 

This  vortex  is  formed  in  the  complicated  flow  along  the  inner  wall. 

The  net  vorticity  in  the  vortex  core  is  a combination  of  the  primary  flow 
vorticity  and  the  secondary  flow  vorticity.  The  primary  vorticity  is 
related  to  the  gradient  of  circulation  along  the  blade  span.  This  vorti- 
city can  be  calculated  by  applying  standard  axisymmetric  flow  theories. 

However,  the  secondary  vorticity  is  due  to  blade  to  blade  pressure 
gradients  and  turning  of  the  flow  with  entering  vorticity  and  can  be 
calculated  from  a set  of  vorticity  transport  equations  applied  along  a 
streamline. 

Experience  has  shown  that  a cavitation  inception  prediction  of  this 
vortex  is  a very  difficult  problem.  All  rotors  operating  with  a boundarv 
layer  have  a vortex  along  the  inner  wall.  The  appearance  of  this 
cavltating  vortex  varies  from  rotor  to  rotor,  as  shown  in  Figure  1.  The 
critical  cavitation  number  can  vary  as  much  as  an  order  to  magnitude. 
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Small  variations  in  the  inlet  velocity  gradient  to  the  rotor  can  cause  a 
significant  change  in  the  critical  cavitation  number. 

1.2  Previous  investigations 

There  has  been  no  systematic  work  on  the  critical  cavitation 
characteristics  of  the  secondary  flow  generated  vortex.  Gearhart  / i / 
investigated  the  secondary  flow  generated  along  the  inner  wall  of  a 
rotor;  however,  no  attempt  was  made  to  correlate  calculations  with  the 
structure  of  the  vortex  or  with  the  critical  cavitation  number  of  the 
vortex.  Nevertheless,  Gearhart  was  one  of  the  first  to  show  the  impor- 
tance of  secondary  flows  in  a rotor  operating  within  a boundary  layer. 

Work  has  been  published  on  the  subjects  of  vortex  cavitation,  the 
structure  of  a vortex,  and  secondary  flows  in  a rotor. 

The  most  extensive  previous  work  on  the  subject  of  vortex  cavitation 
was  conducted  by  McCormick  /_  2 / at  the  Applied  Research  Laboratory 
employing  various  families  of  finite  hydrofoils  which  produced  tip  vortex 
cavitation.  It  was  found  that  the  critical  cavitation  number  depended 
upon  the  boundary  layer  thickness  on  the  pressure  surface  of  the  wing  tip 
and  that  the  thickness  of  the  vortex  core  is  apparently  not  determined  by 
the  induced  drag  of  the  wing.  For  a given  wing  shape,  the  critical  cavi- 
tation number  increased  with  Reynolds  number  and  angle  of  attack  and  was 
nearly  independent  of  aspect  ratio.  Although  McCormick's  investigation 
indicated  the  importance  of  the  boundary  layer,  the  minimum  pressure 
coefficient  and  the  structure  of  the  vortex  were  not  measured.  Further- 
more, although  some  exploratory  investigations  of  the  apparent  nonvaporous 
cavitation  were  conducted,  the  observed  cavitation  was  not  related  to 
measured  values  of  free  gas  content. 
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More  recently,  Hashimoto  j_  3/  has  Investigated  developed  vortex 
cavitation  produced  by  a set  of  stationary  vanes.  This  Investigation  was 
devoted  primarily  to  the  study  of  the  oscillatory  nature  of  the  cavity 
surface  and  a theory  was  developed  to  predict  this  behavior. 

There  is  a large  volume  of  literature  on  the  subject  of  vortex 
structure.  Vortices  have  been  investigated  in  wakes,  over  the  leading 
edge  of  highly-swept  wings,  in  the  wake  of  any  wing  of  finite  span,  in 
pipe  flows,  and  in  nature  as  hurricanes  and  tornadoes.  A major  result  is 
that  these  vortices  are  all  similar  to  the  Rankine  combined  vortex  model 
having  a rotational  region,  which  is  assumed  to  be  of  finite  size, 
outside  of  which  the  motion  is  irrotational . 

No  general  method  for  calculating  the  complete  vortex  structure 
exists,  partly  because  the  conditions  in  which  they  occur  are  varied  and 
because  the  governing  differential  equations  are  nonlinear.  On  the  otner 
hand,  all  vortices  possess  a roughly  axisymmetric  core  and  the  structure 
of  such  cores  can  be  analyzed.  Most  investigators,  such  as  Hall  ]_  4_7 , 
Gartshore  j_  5 _/,  and  Bossel  j_  6_/  deal  with  the  development  of  a vortex 
core  using  boundary  layer  type  equations  given  the  initial  vortex 
structure. 

A considerable  amount  of  vortex  information  has  been  obtained  from 
experimental  observations.  Sarpkaya  j_  7_/ , Harvey  / 8_/,  and  Kirkpatrick 
J_  9_/  have  found  that  under  certain  flow  conditions  the  vortex  can 
abruptly  change  its  structure.  In  studying  swept  wings,  Peckham  and 
Atkinson  [_  10 _/  observed  that  the  core  of  a leading  edge  vortex  expanded 
at  some  distance  downstream  from  the  apex  of  the  wing.  Other  investi- 
gators, such  as  McCormick  l_  11 / , observed  this  expansion  of  the  vortex 

core  in  the  vortex  of  a straight  wing.  This  phenomenon  of  vortex 


breakdown  results  in  an  abrupt  change  in  the  vorticity  distribution  of 
the  vortex  and  a sudden  expansion  of  the  vortex  core. 
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Breakdown  has  been  observed  in  various  form.,  under  quite  different 
conditions.  There  are  several  factors  which  influence  vortex  breakdown. 
The  important  factors  are  (1)  swirl  angle  defined  as  flic  angle  between 
the  axial  velocity  vector  and  the  maximum  tangential  velocity  vector  and 

(2)  an  adverse  pressure  gradient.  In  1972,  Hall  < 1 2 / published  a survey 

article  on  this  phenomenon. 

The  importance  of  the  secondary  flows  in  t rotor  have  been  discussed 
by  many  investigators,  such  as  Hawthorne  ]_  13  /,  the  staff  of  NASA  Lewis 

Research  Center  / 1 4 / , Smith  J_  15/,  and  Lakshminarayana  and  Horlock 

/ 16/.  Experimental  results  show  that  the  secondary  flow  created  in  the 
blade  passage  can  significantly  effect  the  outlet  angle  of  the  blade  row. 

A secondary  flow  is  formed  when  a flew  with  a spanwise  velocity 
gradient  or  boundary  layer  is  turned  by  a blade  row.  It  is  simply  the 
cross-flow  component  of  the  skewed  boundary  layer  which  forms  in  the  blade 
passage.  This  resultant  cross-flow  contains  vorticity  aligned  in  the 
streamwise  direction. 

There  are  no  simple  techniques  tha'’  accurately  predict  these  inner 
wall  velocity  gradient  effects  for  rotational  flows.  Theoretical  under- 
standing of  secondary  flows  through  a rotor  u aided  by  approximate 
solutions  of  the  fluid  flow  equations  that  govern  the  flow  process.  There 
are  many  different  approaches  to  the  derivation  of  secondary  flow  equa- 
tions; however,  Horlock  ]_  17/,  in  a discussion  of  a paper  on  secondary 
vorticity  in  axial  compressor  blade  rows,  shows  that  the  different 
approaches  all  lead  to  essentially  the  same  result,  A review  of  secondary 
flow  theories  can  bo  found  in  the  papers  of  Lakshminarayana  and  Horlock 
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J_  16,  18  /,  the  staff  of  NASA  Lewis  Research  Center  J_  14 / , Hawthorne  and 

Novak  [_  19  / , and  Salvage  J_  20 _/  . 

1 . 3 General  Scope  of  Problem 

The  vortex  formed  near  the  inner  wall  of  a rotor  operating  within  a 
boundary  layer  is  studied  to  determine  the  major  factors  which  influence 
the  minimum  pressure  coefficient  associated  with  the  vortex.  This  vortex 
exists  in  the  low  pressure  region  near  the  inner  wall  of  the  complicated 
flow  behind  a rotor.  In  this  region,  there  are  many  factors  which  can 
influence  the  minimum  pressure  coefficient.  However,  previous  experi- 
mental results  indicate  that  the  vorticity  that  exists  in  the  blade 
passage  controls  the  behavior  of  the  vortex.  This  distributed  passage 
vorticity  appears  to  organize  the  surrounding  low  momentum  fluid  into  a 
vortex. 

This  investigation  of  vortex  cavitation  is  both  experimental  and 

theoretical.  An  inviscid  secondary  flow  analysis  is  developed  which  Is 

employed  to  assess  the  effect  of  passage  streamwise  vorticity  on  the 

circulation  and  core  size  of  the  vortex.  This  analysis  does  not  give  the 

absolute  value  of  C . However,  a correlation  is  made  between  the 
^min 

passage  vorticity  and  vortex  cavitation  data. 

The  experimental  part  of  this  investigation  is  presented  before  the 
theoretical  development  because  a picture  defines  the  vortex  problem 
better  than  words.  An  understanding  of  the  flow  field  is  necessary  in 
order  to  follow  many  of  the  assumptions  used  in  the  theoretical  develop- 


ment . 


2.  BASIC  CONCEPTS 


2 . 1 Various  Rep  lines  of  Cavitation 

Cavitation  flows  are  delineated  by  the  cavitation  number,  ■),  defined 
as 

P -P 

00  v ,,  , 

0 - , (1) 

l/2pVoo2 

where  P , V , P , and  p are  the  pressure  at  infinity,  the  velocity  at 

QO  OO  V 

infinity,  the  vapor  pressure,  and  the  mass  density  of  the  liquid,  respec- 
tively. The  fluid  properties  correspond  to  the  bulk  temperature  of  the 

liquid . 

The  flow  regime  of  particular  concern  is  limited  cavitation,  i.e., 
the  extent  of  cavitation  is  minimized.  The  problem  of  limited  cavitation 
has  been  reviewed  in  recent  years  by  Holl  / 21/  and  lloll,  Arndt,  and 
Billet  122 J. 

The  particular  value  of  a corresponding  to  limited  cavitation  is  the 
limited  cavitation  number,  o^,  or  critical  cavitation  number  given  by 
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where  a,  is  the  desinent  cavitation  number  / 23  /.  The  limited  cavita- 

d — — 

tion  number  can  also  be  determined  by  an  inception  test,  provided 
hysteresis  effects  are  not  involved  j_  24 / . 

There  are  two  general  types  of  limited  cavitation,  namely,  vaporous 
and  nonvaporous  cavitation.  To  visualize  the  various  types  of  limited 
cavitation,  it  is  useful  to  imagine  that  one  is  observing  the  growth  of  a 
single  bubble  as  the  pressure  is  suddenly  reduced.  Vaporous  cavitation 
occurs  at  pressures  less  than  vapor  pressures  and  is  characterized  by  the 
explosive  growth  of  a bubble  due  to  the  rapid  conversion  of  liquid  to 
vapor  at  the  bubble  wall.  There  are  two  types  of  nonvaporous  cavitation, 
namely,  pseudo  and  gaseous  cavitation.  Pseudo  cavitation  occurs  when  a 
bubble  merely  expands  due  to  a reduction  in  pressure  with  the  mass  of  gas 
in  the  bubble  essentially  remaining  constant.  Gaseous  cavitation  occurs 
when  a bubble  grows  in  an  oversaturated  liquid  due  to  the  transport  of 
gas  across  the  interface.  In  contrast  to  vaporous  cavitation,  nonvaporous 
cavitation  can  occur  at  pressures  greater  than  vapor  pressure. 

Of  particular  importance  in  the  study  of  limited  cavitation  is  the 

minimum  pressure  coefficient,  C , given  by 

pmin 


Pmin  l/2pV  2 

CO 

where  P . is  the  minimum  pressure, 
mm  r 

In  general,  results  indicate  for  vaporous  limited  cavitation  that 
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In  principle,  nonvaporous  cavitation  can  occur  below  or  above  vapor 
pressure  /_  21  / . However,  from  a practical  point  of  view,  it  is  the 
latter  type  which  is  the  most  interesting  because  of  the  large  values  oi 
cavitation  numbers  which  are  often  typical  of  these  cases  / 21  /.  Thus, 
for  nonvaporous  limited  cavitation  occurring  at  pressures  above  vapor 
pressure , 


o 


l 


-C 

^min 


(6) 


In  1960,  Holl  ]_  23_/  employing  an  equilibrium  theory  which  ignored 
surface  tension  effects  predicted  for  nonvaporous  cavitation  that 


a.  - -C  + 
*■  P 


G 


l/2pV 


(7) 


where  C is  a characteristic  pressure  coefficient  and  P is  a pressure 
P G 

typical  of  the  noncondensable  gas  in  the  bubbles.  Holl  indicated  that 
there  was  an  upper  limit  to  the  gas  pressure,  namely,  the  equilibrium 
value  given  by  Henry's  Law  which  may  be  expressed  in  the  form 

PG  = , (8) 

max 

where  a is  the  dissolved  gas  content  and  6 is  the  Henry's  Law  constant. 

Several  general  implications  come  from  Equation  (7) . Firstly,  the 

equation  predicts  that  is  a monotonically  decreasing  function  of 

velocity  for  constant  values  of  P_  and  C . Secondly,  the  equation 

G p 

predicts  excessively  large  values  of  at  low  velocities.  These  results 
suggest  that  if  one  is  to  avoid  confusing  nonvaporous  for  vaporous 
cavitation,  it  is  advisable  to  conduct  tests  at  high  velocities.  The 
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aforementioned  tendency  for  to  decrease  with  velocity  for  nonvaporous 
cavitation  is  characteristic  of  data  reported  by  several  authors,  namely, 

Haminitt  et  al.  _/  25  /,  Lindgren  and  Johnson  / 26 / , McCormick  / 1 1 / , and 

Ripken  and  Killen  l_  27 / . 

Thus,  it  is  very  important  to  understand  the  differences  between 
vaporous  and  nonvaporous  cavitation  in  order  to  interpret  test  data.  It 
is  also  apparent  in  view  of  Equations  (5)  and  (6)  that  it  is  very  impor- 
tant to  know  C in  order  to  differentiate  between  various  flow  states. 

P • 
min 

However,  in  many  flows  C data  are  not  available  so  that  other  means 

p . 
mm 

for  differentiating  between  flow  states  are  needed. 

2 . 2 Minimum  Pressure  Coefficient  of  a Vortex 

It  is  appropriate  to  find  a simple  description  of  a vortex  before 
making  a calculation  of  its  minimum  pressure  coefficient.  Unfortunately, 
vortices  are  composed  of  a finite  number  of  vortex  filaments  and  a dif- 
ficulty arises  in  specifying  this  number.  The  effect  of  streamwise 
filaments  is  to  induce  a swirl  velocity  distribution  (V  ) . As  a result 

O 

of  this  swirl  velocity,  a radial  static  pressure  gradient  developes. 
Therefore,  the  direct  method  for  obtaining  the  minimum  pressure  is  to 
numerically  integrate  the  velocity  distribution  from  infinity  to  the 
vortex  axis.  The  minimum  pressure  occurs  at  the  center  of  a vortex. 

The  equation  used  in  the  calculation  of  the  pressure  in  a vortex  for 
axisymmetric  flow  is  the  radial  momentum  equacion, 
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Kor  this  equation,  the  pressure  is  a function  of  the  components  of 
velocity  in  two  planes.  The  component  of  velocity  which  is  tangent  to 

i 

the  streamline  and  is  projected  onto  the  meridional  plane  is  called  the 
meridional  velocity  (V  ) . This  component  of  total  velocity  is  related  to 
the  axial  component  by  the  cosine  of  the  streamline  angle  (p) . The 
second  component  of  velocity  is  known  as  the  tangential  velocity  (V  ) . 

It  is  located  in  the  circumferential  direction  perpendicular  to  the 
meridional  plane. 

As  shown  in  Equation  (9),  there  are  three  contributing  factors  in 
the  development  of  a pressure  gradient  in  a vortex.  The  first  term  is 
directly  related  to  the  centrifugal  force  which  a rotating  fluid  element 
experiences.  The  second  term  is  due  to  a moving  fluid  element  being 
subjected  to  the  streamline  curvature  in  the  meridional  plane.  The 
magnitude  of  this  terra  depends  greatly  on  the  radius  <>!  curvature  of  the 
streamline.  The  third  term  is  due  to  the  convective  acceleration  as  the 
area  of  the  vortex  either  converges  or  diverges. 

If  the  Rankine  combined  vortex  model  is  employed,  namely,  a simple 
rotational  core  combined  with  an  irrotational  outer  flow,  the  pressure  in 
the  center  can  be  calculated  by  integrating  Equation  (9).  For  this  case, 
the  vortex  is  axisyraraetric  with  an  infinite  radius  of  curvature  in  the 
meridional  plane  and  a small  gradient  of  meridional  velocity.  The  result 
is 
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where,  in  the  rotational  core  (r^O , the  tangential  velocity  is  given 

by  V.  = air  and,  in  the  irrotational  region  (r>r  ),  the  tangential 
0 c 

r 

velocity  is  V^  = In  these  relationships,  a>  is  the  rotational 

velocity,  F is  the  circulation,  and  r^  is  the  core  radius. 

The  circulation  can  also  be  related  to  the  vorticity  in  the  vortex 
as 
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where  is  now  the  vorticity  in  the  direction  of  the  axis  of  the  vortex. 
The  vorticity  in  the  axial  direction  can  be  expressed  as 
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where,  for  axisymmetric  flow,  the  derivatives  with  respect  to  the 
O-direction  are  zero. 


2 . 3 Characteristics  of  Vortex  Cavitation 

Vortex  cavitation  is  one  of  the  least  understood  forms  of  cavitation 
It  is  felt  that  this  lack  of  understanding  is  due  to  an  inadequate  knowl- 
edge of  the  flow  field  and  due  to  a confusion  on  types  of  cavitation, 
i.e.,  vaporous  versus  nonvaporous  cavitation.  The  flow  field  can  be 
drastically  altered  by  vortex  breakdown  and,  thus,  influence  the  cavita- 
tion number  describing  the  observed  state  of  cavitation.  Vortex  flows 
tend  to  be  good  collectors  of  gas  bubbles  which  can  cause  nonvaporous 
cavitation.  This  often  leads  to  confusing  nonvaporous  for  vaporous 
cavitation  giving  high  cavitation  numbers. 


3.  EXPERIMENTAL  INVESTIGATION'  OF  THE  FLOW  FIELD 


3.1  Introduction 

Previous  experimental  results  indicate  that  small  changes  in  the 
inlet  velocity  gradient  near  the  inner  wall  of  the  rotor  created  rather 
large  differences  in  the  cavitation  number  of  the  vortex.  In  addition, 
only  small  changes  in  the  thrust  or  torque  of  the  rotor  were  found  tor 
most  changes  in  inlet  velocity  gradient.  These  experimental  result, 
showing  a dependence  on  inlet  velocity  gradient  led  to  this  investigation 
of  secondary  flows. 

Secondary  flows  generated  in  a rotor  for  incompressible  flow  are 
primarily  a function  of  the  inlet  velocity  gradient.  In  order  to 
maximize  the  scale  of  the  inlet  velocity  gradient  for  measurement,  a 
surface  having  an  unfavorable  pressure  gradient  in  front  of  the  rotor 
plane  was  used  in  the  experimental  program.  The  rotor  was  located  near 
the  end  of  this  surface. 

Thus,  the  inlet  velocity  gradient  to  the  rotor  was  due  in  a large 
part  to  the  boundary  layer  on  a surface.  This  boundary  layer  thickness 
in  the  rotor  plane  without  the  rotor  was  the  same  order  of  magnitude  as 
the  rotor  radius.  Also,  the  boundary  layer  thickness  near  the  end  of  this 
surface  was  the  same  order  of  magnitude  as  the  surface  curvature. 

The  experimental  program  used  several  flow  facilities  in  the 
Garfield  Thomas  Water  Tunnel  Building  of  the  Applied  Research  Laboratory 
at  The  Pennsylvania  State  University.  In  all,  there  were  four  different 
test  series  which  were  needed  to  establish  the  character  of  the  flow  field 
near  the  rotor  plane.  Three  of  the  test  series  were  conducted  in  the 
48- inch  water  tunnel.  An  additional  test  was  conducted  in  the  48-inch 
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wind  tunnel.  A description  of  these  flow  facilities  is  given  in 
Reference  j_  28  /. 

To  minimize  any  effects  produced  by  tunnel  wall  interference,  liners 
were  used  in  the  water  tunnel  test  section  and  in  the  wind  tunnel  test 
section.  The  resulting  inner  contour  in  the  test  section  was  determined 
by  a potential  flow  solution  for  an  expanded  surface.  Therefore,  the 
liner  was  a stream  surface  where  the  surface  is  in  a flow  to  infinite 
extent . 

Each  test  series  will  not  be  discussed  individually  because  in  most 
cases  parts  of  several  tests  were  needed  to  establish  one  data  curve.  An 
error  analysis  of  the  experimental  data  is  given  in  Appendix  A.  The 
experimental  results  of  the  flow  field  measurements  were  divided  into  two 
main  groups.  The  first  group  consists  of  measurements  of  the  primary 
flow  field,  such  as  rotor  inlet  and  outlet  velocity  profiles.  Data 
obtained  to  establish  the  character  of  the  vortex  are  the  second  group. 
This  includes  vortex  velocity  measurements,  pressure  measurements  along 
the  surface,  and  vortex  cavitation  data. 

3. 2 Basic  Flow  Configurations 

There  were  several  ways  to  vary  the  inlet  velocity  gradient  to  the 
rotor.  The  easiest  one  was  to  change  the  free  stream  velocity  (V  ) for  a 
specified  rotor  flow  coefficient  defined  as 

+ = TT1-  , (13) 

TIP 

where  U p is  the  rotor  tip  velocity.  This  gives  the  influence  of 
Reynolds  number  on  the  secondary  flow  generated  vortex. 
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More  important  than  the  aforementioned  cases  are  inlet  velocity 
gradient  variations  for  the  same  Reynolds  number  based  on  free  stream 
velocity  and  rotor  diameter.  These  variations  were  accomplished  by 
adding  a screen  to  the  surface  upstream  of  the  rotor,  operating  the  rotor 
at  different  flow  coefficients,  and  adding  struts  to  the  surface.  Adding 
a screen  to  the  surface  increased  the  thickness  of  the  boundary  layer 
entering  the  rotor.  The  upstream  struts  consisted  >f  four  struts  placed 
at  the  0°,  90°,  180°,  and  270°  points  on  the  surface  which  made  the  inlet 
velocity  profile  nonaxisymmetric . 

Each  variation  in  inlet  velocity  gr.  Jit-nt  at  a constant  Reynolds 
number  was  a basic  flow  configuration.  A list  of  the  basic  flow  configu- 
rations is  given  in  Table  1.  In  many  cases,  the  free  stream  velocity  was 
also  varied  for  a basic  flow  configuration.  These  thirteen  basic  flow 
configurations  formed  the  test  matrix  for  the  experimental  program. 

3 . 3 Primary  Flow  Field  Measurements 

3.3.1  Surface  Velocity  Profile  Measurements.  Several  velocity 
profiles  were  measured  in  the  plane  of  the  rotor  without  the  rotor 
installed.  In  the  wind  tunnel  at  a Reynolds  number  of  3.0  x 105  based  on 
rotor  diameter  and  free  stream  velocity,  measurements  were  obtained  for 
(1)  Basic  Flow  No.  11  (without  upstream  struts,  without  upstream  screen, 
and  without  rotor),  (2)  Basic  Flow  No.  12  (without  upstream  struts,  with 
upstream  screen,  and  without  rotor),  and  (3)  Basic  Flow  No.  13  (with 
upstream  struts,  without  upstream  screen,  and  without  rotcr) . In  the 
water  tunnel  at  a Reynolds  number  of  8.8  x 10s,  measurements  were  again 
obtained  for  Basic  Flow  No.  11.  Only  Basic  Flow  No.  11  is  amenable  to 
theoretical  boundary  layer  analysis. 
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TABLE  1 

BASIC  FLOW  CONFIGURATIONS 


Basic 

Flow  Nos. 

Upstream 

Struts** 

Upstream 

Screen 

Flow 

Coefficient 

<f>*** 

Rotor 

Cap 

1 

Without 

Without 

Design 

Conical 

2 

Without 

Without 

0.9  Design 

Conical 

3 

Without 

With 

Design 

Conical 

4 

With 

Without 

Design 

Conical 

5 

With 

Without 

1.1  Design 

Conical 

6 

With 

Without 

0.9  Design 

Conical 

7 

With 

With 

Design 

Conical 

8 

With 

With 

0.9  Design 

Conical 

9 

Without 

With 

0.9  Design 

Conical 

10 

With 

Without 

Design 

Truncated 

11 

Without 

Without 

No  Rotor 

Conical 

12 

Without 

With 

No  Rotor 

Conical 

13 

With 

Without 

No  Rotor 

Conical 

14* 

Without 

Without 

Design  at  Rerej 

Conical 

15* 

Without 

Without 

Design  at  3 Re  ^ 

Conical 

16* 

Without 

Without 

0.9  Design  at  Re  . 

ref 

Conical 

17* 

Without 

Without 

0.9  Design  at  3 Re  c 
ref 

Conical 

These  are  special  cases  for  theoretical  calculations  only. 

**The  upstream  struts  consisted  of  four  struts  placed  at  0°,  90°,  180°, 
and  270  on  the  surface. 

Vo 

***The  flow  coefficient  is  defined  as  = — . 

TIP 
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For  the  wind  tunnel  measurements,  the  local  velocity  was  calculated 
from  total  pressure  measurements  after  the  static  pressure  distribution 
through  the  boundary  layer  was  established.  The  accuracy  of  these 
measurements  is  discussed  in  Appendix  A.  In  all,  three  total  pressure 
rakes  and  one  static  pressure  rake  were  used.  As  the  rakes  were  rotated, 
pressures  were  recorded  at  2°  increments  for  a total  of  360  In 
addition,  static  pressures  were  obtained  from  pressure  taps  located  on 
the  suriace  to  insure  the  true  static  pressure  near  the  inner  wall. 
Further  details  of  these  measurements  are  given  in  Reference  j_  29/. 

The  velocity  profiles  are  shown  in  Figure  2 for  Basic  Flow  No.  11, 
in  Figure  3 for  Basic  Flow  No.  12,  and  in  Figure  4 for  Basic  Flow  No.  13. 
The  mean  velocity  profiles  shown  in  the  figures  were  obtained  by  curve 
fitting  the  data  points.  Also,  the  data  points  represent  the  circumfer- 
ential average  of  180  experimental  points.  Tables  listing  the  mean 
velocity  profiles  are  given  on  pages  26,  27,  and  28  of  Reference  j_  29_/. 

The  effect  of  adding  an  upstream  screen  can  be  seen  by  comparing 
Figure  2 to  Figure  3.  The  screen  essentially  thickens  the  momentum 
deficit  near  the  inner  wall.  The  effect  of  adding  upstream  struts  can  be 
seen  by  comparing  Figure  2 to  Figure  4.  This  comparison  shows  that  the 
upstream  struts  add  momentum  to  the  deficit  near  the  inner  wall  and 
increase  the  local  velocity.  This  resulting  deficit  is  controlled  by  the 
secondary  flows  created  at  the  intersection  of  the  struts  and  the  inner 
wall.  A more  detailed  description  of  the  circumferential  velocity  varia- 
tion due  to  the  struts  is  given  in  the  next  section. 

3.3.2  Rotor  Inlet  Velocity  Measurements.  Several  inlet  velocity 
profiles  were  measured  approximately  one-half  chord  length  in  front  of 


the  rotor.  These  measurements  were  done  in  the  48-inch  wind  tunnel  at  a 


VELOCITY  RATIO,  V/V 

m co 

Figure  2,  Meridional  Velocity  Profile  for  Basic  Flow  No,  11 
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Figure  4.  Meridional  Velocity  Profile  for  Basic  Flow  No.  13 
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Reynolds  number  of  3.0  x 10J  using  several  pressure  rakes  and  one  static 
pressure  rake.  A schematic  showing  the  plane  of  the  measurements  and  the 
location  of  the  rotor  is  shown  in  Figure  5.  The  accuracy  of  these 
measurements  is  discussed  in  Appendix  A. 

The  mean  incoming  velocity  profiles  obtained  by  curve  fitting  the 
data  points  are  shown  in  Figures  6 through  10.  Also,  each  data  point 
represents  the  circumferential  average  of  180  experimental  points.  The 
following  is  the  list  of  flow  configurations  for  which  the  inlet  velocity 
profiles  were  obtained: 

Figure  6 - Basic  Flow  No.  1 

Figure  7 - Basic  Flow  No.  2 

Figure  8 - Basic  Flow  No.  3 

Figure  9 - Basic  Flow  No.  4 

Figure  10  - Basic  Flow  No.  7 

A comparison  of  the  circumferentially  averaged  velocity  measurements 
obtained  with  the  rotor  (Figure  6)  and  without  the  rotor  (Figure  2)  show, 
that  the  mean  velocities  are  higher  for  the  case  with  the  rotor.  This  is 
particularly  true  near  the  inner  wall.  It  is  probable  that  the  higher 
mean  velocities  found  with  the  rotor  are  caused  by  streamline  convergence. 
This  effect  is  produced  by  the  favorable  pressure  gradient  generated  by 
the  rotor. 

The  effect  of  an  upstream  screen  on  the  rotor  inlet  velocity  profile 
can  be  seen  by  comparing  Figure  6 to  Figure  8.  The  screen  essentially 
thickens  the  momentum  deficit  near  the  inner  wall  of  the  rotor  inflow 
region. 

The  effect  of  upstream  struts  can  be  seen  by  comparing  Figure  6 to 
Figure  7.  The  results  are  similar  to  those  previously  discussed  for 
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Figure  6.  Rotor  Inlet  Velocity  Profile  for  Basic  Flow  No.  1 
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Figure  7.  Rotor  Inlet  Velocity  Profile  for  Basic  Flow  No.  2 


NONDIMENSIONAL  DISTANCE  FROM  SURFACE, 


d(V  /V  > 
m 


NONDIMENSIONAL  NORMAL  VORTICITY.  un  = d(y/R 


O.S 


1.0 


1.5 


2.0 


2.5 


R 

3.0 


3.5 


Figure  9. 


Rotor  Inlet  Velocity  Profile  for  Basic  Flow  No.  4 
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velocity  measurements  obtained  without  the  rotor  and  with/without 
upstream  struts.  However,  the  mean  velocities  are  higher  for  the  cases 
with  the  rotor. 

More  details  of  the  effect  of  the  struts  are  shown  in  Figures  11 
through  16.  These  figures  show  the  circumferential  variation  of  velocity 
at  various  radial  positions.  Even  at  relatively  large  distances  from  the 
surface,  sharp  depressions  in  the  velocity  profile  are  produced.  Near 
the  surface  and  strut  intersection,  vortices  are  formed  which  entrain 
fresh  fluid  into  the  immediate  neighoorhood  of  the  surface  and  displace 
sluggish  fluid  from  the  boundary  layer.  As  a result,  large  peaks  and 
valleys  are  created  in  the  circumferential  variation  of  velocity. 

For  the  many  flow  configurations  considered  in  this  investigation, 
the  inlet  velocity  profile  to  the  rotor  is  basically  a surface  boundary 
layer  profile  with  a pertubation  imposed  by  the  rotor.  Therefore, 

* 

standard  boundary  layer  parameters  such  as  displacement  thickness  (i  ) 
and  momentum  thickness  (0)  might  be  meaningful  parameters  for  describing 
the  inlet  velocity  profile.  However,  the  exact  boundary  layer  thickness 

(5)  defined  as  the  distance  from  the  surface  when  V is  99%  of  V cannot 

x «> 

be  accurately  determined  for  the  cases  witli  the  rotor.  This  is  due  to 
the  influence  of  the  rotor  on  the  inlet  velocity  profile.  Thus,  the  data 
were  analyzed  to  obtain  the  boundary  layer  parameters.  The  results  are 
shown  In  Table  2 using  planar  two-dimensional  definitions  and  in  Table  3 
using  axisymmetric.  definitions. 

The  two-dimensional  results  (Table  2)  were  obtained  by  integrating 
the  following  relationships 
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Figure  11.  Comparison  of  Rotor  Inlet  Velocity  Measured  0.062 
Inch  from  Surface  with  and  without  Upstream  Struts 
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Figure  14.  Comparison  of  Rotor  Inlet  Velocity  Measured  0.439 
Inch  from  Surface  with  and  without  Upstream  Struts 
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Basic 
Flow  No. 

Upstream 

Struts 

Upstream 
Screen 

Flow 

Coefficient 

6 

I nches 

6 

1 nches 

H 

I ! 

Without 

Without 

Design 

0 . 644 

1.179 

1.39 

1 

Without 

Without 

0.9  Design 

1.414 

1.057 

1.34 

I 3 

Without 

With 

Design 

2.003 

1.495 

1.34 

4 

With 

Wi  thout 

Design 

1.825 

1.328 

1.37 

7 

With 

With 

Design 

2.140 

1.574 

1.36 

TABLE  3 

AX I SYMMETRIC  BOUNDARY  LAYER  ANALYSIS 


6 - displacement  thickness 

0 - momentum  thickness 

6 - boundary  layer  thickness 

k 

H - = A /0 
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where  Y is  the  normal  distance  from  the  inner  wall. 

On  the  other  hand,  the  incoming  velocity  gradient  is  axisymmetric  so 
the  appropriate  axisymmetric  definitions  for  displacement  and  momentum 
thickness  are 
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Patel  j_  30 _/  has  shown  that  there  is  a relationship  between  planar  defini- 
tions and  the  axisymmetric  definitions  so  that  the  calculation  of  only  one 
set  of  parameters  is  necessary.  The  calculations  were  done  using  the  two- 
dimensional  relationships. 

The  results  of  using  boundary  layer  definitions  parallel  the 
conclusions  from  the  figures,  i.e.,  (1)  the  screen  increases  the  momentum, 
thickness,  (2)  the  average  effect  of  the  upstream  struts  increases  the 
momentum  thickness,  and  (3)  operating  the  rotor  at  less  than  design  flow 
coefficient  decreases  the  momentum  thickness.  In  all  configurations,  the 
rotor  operated  within  the  estimated  inlet  boundary  layer. 

The  slope  of  the  inlet  velocity  profile  to  the  rotor  determines  the 
amount  of  vorticity  entering  the  rotor  for  the  axisymmetric  case.  The 
nondimensional  normal  vorticity  is  defined  as 


0) 
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d(v  /v  ) 

m 00 
d(Y/RR) 


(16) 
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where  V /V  is  a nondimens ional  meridional  velocity  profile,  Y is  the 
m •»' 

normal  distance  from  the  inner  wall,  and  R is  the  radius  of  the  rotor. 

R 

The  radius  of  the  rotor  was  chosen  as  the  normal i zing  parameter  because 
of  the  uncertainty  of  the  boundary  layer  thickness. 

The  normal  vorticity  data  obtained  from  the  velocity  profiles  art- 
shown  in  Figures  6 through  10  and  are  listed  in  Tables  4 through  8 in 
Reference  / 29/.  The  important  conclusion  is  that  for  the  average  inlet 
profile  to  the  rotor  the  upstream  struts  reduce  the  amount  of  normal 
vorticity  entering  the  rotor  as  compared  to  the  no  strut  cases. 

3.3.3  Rotor  Outlet  Velocity  Measurements . The  velocity  profiles 
behind  the  rotor  were  measured  for  several  rotor  operating  flow  coef- 
ficients and  several  upstream  surface  configurations.  The  location  of 
these  measurements  was  approximately  one-half  chord  length  downstream  of 
the  rotor  exit  plane  as  shown  in  Figure  5.  The  results  were  obtained  at 
a Reynolds  number  of  8.8  x 10-  based  on  rotor  diameter  in  the  ARL  48-inch 
diameter  water  tunnel.  Most  of  the  velocity  component  data  were  obtained 
by  using  a 5-hole  prism  probe  (see  Reference  / 31  /);  however,  a Laser 
Doppler  Anemometer  was  used  to  verify  some  of  the  rotor  outlet  5-hole 
probe  measurements.  The  calibration  of  the  5-hole  prism  probe  used  in 
the  measurements  and  a more  detailed  description  of  the  experimental 
program  is  given  in  Reference  j_  32 _/.  In  addition,  the  accuracy  of  these 
measurements  is  discussed  in  Appendix  A. 

The  mean  velocity  profiles  are  shown  in  Figures  17  through  25  and 
are  listed  in  Tables  2 through  10  in  Reference  / 32  /.  The  following  is 
the  list  of  flow  configurations  for  which  profiles  were  obtained: 

Figure  17  - Basic  Flow  No.  1 


Figure  18  - Basic  Flow  No.  2 
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Figure  19  - Basic  Flow  No.  3 
Figure  20  - Basic  Flow  No.  4 
Figure  21  - Basic  Flow  No.  5 
Figure  22  - Basic  Flow  No.  6 
Figure  23  - Basic  Flow  No.  7 
Figure  24  - Basic  Flow  No.  8 
Figure  25  - Basic  Flow  No.  9 
The  Basic  Flow  Nos.  are  defined  in  Table  1. 

These  figures  show  both  data  obtained  using  a 5-hole  prism  probe  and 
a LDA.  The  LDA  was  used  primarily  to  determine  how  close  measurements  of 
the  tangential  velocity  (V^)  could  be  made  to  the  rotating  surface  before 
measuring  the  effect  of  rotation  on  the  flow.  Measurements  were  made 
without  the  rotor  but  with  the  rotating  cap.  The  tangential  velocity 
ratio  (Vg/V^)  was  found  to  increase  dramatically  to  a value  of  0.61  for 
design  flow  coefficient  within  a distance  of  0.20  inch  from  the  surface. 
Measurements  of  the  tangential  velocity  profile  with  the  rotor  were  not 
made  within  a distance  of  0.20  inch  from  the  surface.  The  tangential 
velocity  profile  near  the  inner  wall  is  a result  of  the  secondary  flows 
and  rotor  loading.  In  all  cases,  the  normalized  axial  velocity  ratio 
(V^/V^)  has  a parabolic  profile  downstream  of  the  rotor  as  a result  of  the 
loading  of  the  rotor.  This  outlet  axial  profile  is  compared  to  the 
profile  that  would  exist  without  the  rotor  and  the  inlet  velocity  profile 
to  the  rotor  in  Figure  26  for  Basic  Flow  No.  1.  In  this  figure,  the 
velocity  profile  without  the  rotor  and  the  inlet  velocity  profile  were 
corrected  from  the  meridional  velocity  to  the  axial  velocity.  The  change 
in  velocity  was  small  because  the  streamline  angles  were  less  than  8 


degrees . 


Figure  19.  Rotor  Outlet  Velocity  Profiles  for  Basic  Flow  No.  3 
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Figure  20.  Rotor  Outlet  Velocity  Profiles  for  Basic  Flow  No.  4 
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Figure  23.  Rotor  Outlet  Velocity  Profiles  for  Basic  Flow  No.  7 
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Figure  24.  Rotor  Outlet  Velocity  Profiles  for  Basic  Flow  No.  8 
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Figure  25.  Rotor  Outlet  Velocity  Profiles  for  Basic  Flow  No.  9 
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Figure  26.  Comparison  of  Velocity  Profiles  for  Basic  Flow  No.  1 


50 


A comparison  of  tangential  velocity  ratios  (V^/V)  obtained  with  the 

rotor  operating  at  10%  high  in  design  flow  coefficient  (Figure  21),  design 

flow  coefficient  (Figure  20),  and  10%  low  in  design  flow  coefficient 

(Figure  22)  shows  that,  as  the  flow  coefficient  increases,  the  amount  of 

turning  decreases.  Also,  the  axial  velocity  ratio  (V  (V  ) decreases  with 

x “ 

increasing  flow  coefficient.  This  is  particularly  true  near  the  inner 
wall . 

The  influence  of  upstream  struts  on  the  measured  profiles  can  be 
noted  by  comparing  Figure  17  with  Figure  20  and  Figure  22  with  Figure  IS. 
It  shows  that  adding  upstream  struts  causes  an  increase  in  both  V„/V  and 

9 cr. 

V /V  near  the  inner  wall  of  the  rotor. 

X 00 

3.4  Vortex  Flow  Field  Measurements 

3.4.1  Vortex  Profile  Measurements  on  Rotor  Cap.  So  far,  only  the 
flow  measurements  near  the  rotor  plane  have  been  presented.  As  can  be 
noted  from  the  measurements  obtained  near  the  rotor  exit  plane  (Figures  1/ 
through  25),  the  tangential  velocity  profile  does  not  show  a well  defined 
vortex  near  the  rotor  inner  wall  in  this  plane. 

As  shown  in  Figure  5,  downstream  of  the  rotor  is  a cap.  A conical 
cap  was  used  because  it  does  not  cause  flow  separation  and  therefore  the 
flow  can  be  calculated  along  the  conical  boundary. 

The  cavitation  associated  with  the  vortex  occurred  near  the  end  of 
this  cap.  Therefore,  flow  measurements  were  obtained  at  the  end  of  the 
cap.  In  addition,  static  pressure  taps  were  located  along  the  cap  and 
these  results  will  be  discussed  in  the  next  section. 

The  behavior  of  the  vortex  when  subjected  to  basic  flow  changes  can 
be  seen  from  measurements  made  at  the  end  of  the  cap.  These  measurements 
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were  made  using  a Laser  Doppler  Anemometer  (LDA) . Preliminary  flow  field 
observations  using  dye  showed  that  a probe  could  not  be  used  to  measure 
the  vortex  structure  near  the  end  of  the  cap.  Introduction  of  any  probe 
into  the  vortex  core  caused  vortex  breakdown  which  resulted  in  the  vortex 
core  diffusing  rapidly  in  front  of  the  probe.  This  led  to  the  development 
of  the  LDA  for  the  48-inch  water  tunnel. 

A picture  of  the  LDA  setup  for  the  43-inch  water  tunnel  is  shown  in 
Figure  27.  The  dual  beam-backscatter  mode  was  used  for  the  measurements 
because  optics  could  not  be  mounted  on  both  sides  of  the  water  tunnel. 

In  the  dual  beam  mode,  the  laser  beam  is  separated  into  two  parallel 
beams  of  equal  intensity.  These  two  beams  cross  and  focus  at  the  focal 
point  of  a lens.  The  scattered  light  originating  from  the  focal  point  is 
collected  through  the  focusing  lens  and  is  directed  to  a photomultiplier 
tube.  A Spectra  Physics  Model  164  Argon-Ion  laser  was  used  as  the 
coherent  light  source.  A Thermo-Systems  Model  1090  Tracker  was  used  to 
process  the  signal. 

The  measurements  of  the  flow  structure  at  the  end  of  the  cap  for 
several  basic  flows  are  shown  in  Figures  28  through  34.  A discussion  of 
the  accuracy  of  these  measurements  is  given  in  Appendix  A.  These  figures 
show  the  complexity  of  the  vortex.  First  of  all,  the  vortex  is  not  steady 
so  that  the  measurements  represent  the  mean  value  at  a given  radial 
position.  In  addition,  the  laser  light  probe  volume  is  of  finite  size 
which  will  also  contribute  to  the  averaging  process.  Secondly,  there  were 
extremely  large  variations  in  the  tangential  velocity  for  some  basic  flows 
which  are  not  shown  in  the  figures.  Variations  as  large  as  twice  the  mean 
tangential  velocity  were  measured. 


Figure  27  - Fixture  of  I. DA  for  48- 1 noli  Water  Tunne 
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Figure  29.  Velocity  Profiles  on  Rotor  Cap  for  Basic  Flow  No.  2 
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Figure  32.  Velocity  Profiles  on  Rotor  Cap  for  Basic  Flow  No.  5 


60 


The  axial  velocity  profile  at  the  end  of  the  cap  shown  in  Figures  28 
through  34  has  a similar  shape  as  the  axial  velocity  profile  measured 
near  the  rotor  exit  plane  as  shown  in  Figures  17  through  25.  However, 
the  tangential  velocity  profile  at  the  end  of  the  rotor  cap  shows  a well 
defined  vortex  which  measurements  near  the  rotor  exit  plane  do  not  show. 

3.4.2  Surface  Pressure  Measurements  Near  the  Rotor  Plane.  Static 
pressures  were  measured  along  the  surface  including  upstream  of  the  rotor 
plane  and  downstream  of  the  rotor  along  the  rotor  cap.  In  all,  five 
pressure  taps  were  in  front  of  the  rotor,  two  were  in  the  blade  passage 
of  the  rotor,  and  five  were  on  the  rotor  cap.  The  axial  location  of  the 
pressure  taps  forward  of  the  end  of  the  rotor  cap  is  given  in  Table  4. 

The  first  seven  pressure  taps  from  the  end  of  the  cap  were  located 
on  a rotating  surface.  These  pressure  measurements  were  obtained  by 
using  a rotating  scani-valve  mounted  at  the  end  of  the  drive  shaft  on  the 
center  line.  This  device  as  shown  in  Figure  35  consists  of  24  pressure 
ports  which  can  be  internally  connected  to  the  outlet  pressure  port.  In 
turn,  this  outlet  pressure  port  rotates  within  a stationary  pressure 
cavity.  This  cavity  is  finally  connected  to  a pressure  transducer. 

Because  the  cap  pressure  taps  are  located  off  the  axis  of  rotation, 
the  measured  pressures  in  the  rotating  frame  will  include  an  additional 
pressure.  Therefore,  the  pressure  measured  on  the  axis  of  rotation  is 

= PT  - l/2p(rT<2)2  , (17) 

where  i is  the  angular  velocity  and  r is  the  distance  from  the  axis  of 
rotation  to  the  tap  location.  The  static  pressure  coefficient  for  the 
surface  pressure  tap  becomes 
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TABLE  4 

LOCATION  OF  SURFACE  PRESSURE  TAPS 


'ap  No. 


Distance  Forward  of  End  of  Rotor  Cap 
( inches) 


On  Stationary  Surface 

P„  15.50 

T 

1 

PT  9.50 

2 

PT  7.90 

3 

PT  5.03 

4 

PT  3.60 

5 

On  Rotor 

P 3.34 

T 

6 

PT  2.71 

7 


On  Rotor  Cap 

PT  0.66 

8 

PT  0.50 

9 

PT  0.34 

10 

PT  0.17 

1 1 

PT  0.00 

12 


63 


r 

i 


i 


! 


VPs 


l/2pV  2 1/20V  2 


(18) 


where  $ is  the  rotor  flow  coefficient  and  R is  the  rotor  radius. 

R 

The  static  pressure  coefficients  measured  on  the  body  surface 
without  the  rotor  are  shown  in  Figure  36.  The  accuracy  of  these  measure- 
ments is  discussed  in  Appendix  A.  The  effect  of  a rotating  cap  on  the 
flow  field  can  be  seen  by  the  decrease  in  static  pressure  coefficient 
near  the  end  of  the  rotor  cap. 

The  effect  of  adding  the  rotor  to  the  surface  is  shown  in  Figure  37. 
In  this  figure,  the  pressure  coefficients  along  the  surface  are  compared 
for  Basic  Flow  No.  11  without  a rotor  but  with  cap  and  Basic  Flow  No.  1 

with  a rotor.  The  rotor  causes  a pressure  rise  at  the  exit  plane  (P^  ), 

7 

and  this  pressure  progressively  decreases  downstream.  An  apparent 
pressure  minimum  occurs  at  the  end  of  the  cap.  In  addition,  the  pressure 
coefficients  upstream  of  the  rotor  plane  are  decreased.  This  indicates 
an  increase  in  velocity  to  the  rotor.  The  influence  of  the  rotor  on  the 
surface  pressure  coefficient  becomes  small  approximately  two  rotor 
diameters  upstream. 

A more  detailed  picture  of  the  pressure  coefficients  on  the  cap 
downstream  of  the  rotor  is  given  in  Figures  38  and  39.  In  all  cases,  the 
pressure  coefficient  decreased  along  the  cap  with  the  measured  minimum  at 
the  end  of  the  cap  and  the  pressure  coefficient  decreased  at  each  tap 
location  as  the  velocity  increased.  The  most  interesting  result  is  the 
effect  of  adding  upstream  struts  (Figure  39)  on  the  pressure  coefficients 
along  the  cap.  In  comparing  Figure  38  to  Figure  39,  the  pressure 
coefficients  are  lower  along  the  cap  without  upstream  struts  with  a 


Static  Pressure  Coefficients  Along  Surface  without  Rotor 


WITHOUT  UPSTREAM  STRUTS 
WITHOUT  SCREEN 

o WITHOUT  ROTOR,  CAP  NOT  ROTATING 
a WITH  ROTOR  (BASIC  FLOW  No.  1) 


DISTANCE  FORWARD  OF  CAP  TIP  (inches) 


PRESSURE  COEFFICIENT,  C 


DISTANCE  FORWARD  OF  END  OF  ROTOR  CAP  (Inches) 

Figure  39.  Static  Pressure  Coefficients  on  Rotor  Cap  for 
Basic  Flow  No.  4 


notable  exception  near  the  end  of  the  rotor  cap.  Near  the  end,  the 
pressure  coefficient  is  lower  with  added  upstream  struts. 

3.4.3  Vortex  Cavitation  Data.  Desinent  cavitation  number  data  were 
obtained  for  different  inlet  velocity  profiles  to  the  rotor,  i.e.,  for 
the  various  Basic  Flows.  In  all  cases,  cavitation  in  the  vortex  system 
occurred  on  the  end  of  the  cap  located  downstream  of  the  rotor.  The 
occurrence  of  the  cavitation  was  sporadic. 

The  air  content  was  approximately  3.1  ppm  for  all  of  the  experiments 
except  for  one  experiment  which  had  a range  of  air  contents.  The  air 
content  of  3.1  ppm  was  chosen  because  gas  effects  are  reduced  and  the 
relative  saturation  level  was  always  much  less  than  unity.  The  abbrevia- 
tion ppm  denotes  moles  of  air  per  million  moles  of  water.  A more 
detailed  description  of  the  cavitation  experiments  is  given  in 
Reference  ]_  33_/. 

The  desinent  cavitation  number  data  are  shown  in  Figures  40  through 
46  where  the  effects  of  the  following  factors  are  displayed: 

Figure  40  - Effect  of  Air  Content 

Figure  41  - Effect  of  Cap  Design 

Figure  42  - Effect  of  Flow  Coefficient 

Figure  43  - Effect  of  Upstream  Struts 

Figure  44  - Effect  of  Upstream  Struts  and  Screens 

Figure  45  - Effect  of  Screen  and  Flow  Coefficient  with  Struts 

Figure  46  - Effect  of  Screens  and  Flow  Coefficient  without  Struts 

Thus,  Figures  40  through  46  show  the  effects  on  the  desinent  cavitation 
number  due  to  air  content  level  and  cap  design  in  addition  to  basic  flow 
conf igurations . 
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WITHOUT  UPSTREAM  STRUTS 

DESIGN  FLOW  COEFFICIENT 
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Figure  46.  Cavitation  Data  - Effect  of  Screen  and  Flow  Coefficient 
without  Struts  (Basic  Flows  Nos.  1,  2,  3,  and  9) 
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In  general,  an  increase  in  gas  content,  the  addition  of  upstream 
struts,  or  a 101  decrease  in  flow  coefficient  causes  the  cavitation 
number  to  increase  for  a constant  Reynolds  number.  In  contrast  to  these 
results,  the  addition  of  an  upstream  screen  or  a 10%  increase  in  flow 
coefficient  causes  the  cavitation  number  to  decrease  for  a constant 
Reynolds  number. 

Data  in  Figure  42  shows  that  a decrease  in  the  flow  coefficient  by 
10%  causes  a dramatic  increase  in  the  cavitation  number,  whereas  a 10% 
increase  in  the  flow  coefficient  causes  the  opposite  trend.  The  cavita- 
tion number  at  the  design  flow  coefficient  increases  slowly  with  velocity, 
whereas  at  10%  below  the  design  flow  coefficient,  it  increases  rapidly 
with  velocity.  In  contrast  to  these  trends,  at  a flow  coefficient  10% 
above  the  design  value,  the  cavitation  number  decreases  with  velocity. 

The  effect  of  upstream  struts  is  shown  by  the  data  in  Figures  43  and 
44,  where  the  struts  are  in  the  form  of  four  struts  upstream  of  the  rotor. 

It  is  seen  that  the  addition  of  the  struts  increases  the  cavitation  number 
by  50  to  70%. 

The  effect  of  cap  design  is  shown  by  the  data  in  Figure  41.  The 

average  of  the  data  indicates  that  the  conical  cap  causes  cavitation 

* 

numbers  wh.iqjh  increase  with  velocity  and  range  between  4.5  and  5.1.  In 
contrast  to? this  trend  a truncated  cap  causes  much  higher  cavitation 

I 

numbers  which  tend  to  decrease  with  an  increase  in  velocity.  This 

i 

velocity  trend  is  in  agreement  with  the  equilibrium  theory  discussed  in 

Chapter  2.  Assuming  that  C and  P are  essentially  constant,  then 

Pmin  G 


Equation  (7)  can  be  expressed  as 
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The  data  for  the  truncated  cap  approximate  this  trend.  Apparently,  the 
separated  flow  increases  the  residence  time  of  bubbles  in  the  low 
pressure  region  which  accentuates  the  gas  effects  predicted  by  the 
equilibrium  theory. 

The  data  in  Figure  40  for  a conical  cap  also  show  gas  effects  on  the 
cavitation  number.  As  the  gas  content  increases,  the  cavitation  number 
tends  to  increase. 

As  discussed  in  Chapter  2,  previous  experimental  results  indicate 
for  vaporous  cavitation  that 


a , = - C , (20) 

d p . 

min 

assuming  that  air  content  effects  are  small.  Therefore,  a comparison  can 

be  made  between  the  measured  pressure  coefficients  at  the  end  of  the  cap 

(C  ) and  the  average  cavitation  data  (a.), 
p , d 

min 

Comparisons  are  given  in  Figure  47  for  basic  flows  without  struts  and 
in  Figure  48  for  basic  flows  with  struts.  A closer  agreement  is  found  for 
the  cases  without  struts  than  with  struts.  However,  the  correlation 
between  the  measured  pressure  coefficients  at  the  end  of  the  rotor  cap  and 
average  cavitation  data  is  at  best  poor.  As  the  free  stream  velocity  is 
increased,  the  correlation  becomes  increasingly  worse. 

3.5  Summary  of  Experimental  Results 

Most  of  this  chapter  is  comprised  of  figures  showing  the  flow  field 
measurements  and  vortex  cavitation  data  for  the  many  basic  flows.  These 
measurements  define  the  flow  field  from  the  rotor  inlet  to  the  end  of  the 
rotor  cap.  The  following  is  a summary  of  measurements  that  were  made: 
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Figure  47.  Comparison  Between  Static  Pressure  Coefficients 
at  End  of  Rotor  Cap  with  Cavitation  Data  for 
Several  Basic  Flows  without  Struts 
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48-Inch  Wind  Tunnel 

(1)  Rotor  inlet  velocity  profiles  were  measured  with  separate 
total  and  static  probes. 

(2)  Static  pressures  were  measured  on  the  stationary  surface 
upstream  of  the  rotor  plane. 

48-Inch  Water  Tunnel 

(1)  Rotor  outlet  velocity  profiles  were  measured  with  a 
5-hole  prism  probe. 

(2)  Static  pressures  were  measured  on  the  stationary  surface 
upstream  of  rotor  plane  and  on  the  rotating  cap  downstream 
of  rotor  plane. 

(3)  Desinent  cavitation  numbers  for  the  vortex  were  obtained 
for  many  basic  flow  configurations. 

(4)  Velocity  profiles  were  measured  at  the  end  of  the  rotor 
cap  using  a LDA. 

(5)  Several  rotor  inlet  and  outlet  velocity  profiles  were 
measured  using  a LDA  as  a check  on  previous  probe 
measurements . 

A discussion  of  the  accuracy  of  these  measurement  instruments  is  given  in 
Appendix  A. 

Desinent  cavitation  numbers  for  the  vortex  were  found  to  vary 
significantly  as  the  inlet  velocity  profile  to  the  rotor  was  changed. 
Changes  to  the  inlet  velocity  profile  were  accomplished  by  (1)  adding  a 
screen  to  the  upstream  surface,  (2)  adding  upstream  struts  to  the  surface, 

(3)  operating  the  rotor  at  different  flow  coefficients,  and  (4)  varying 
the  free  stream  velocity  for  constant  flow  coefficient.  Each  variation  in 
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inlet  velocity  profile  at  a constant  Reynolds  number  was  given  a Basic 
Flow  Number. 

An  analysis  of  the  inlet  velocity  profile  data  shown  in  Figures  6 
through  10  shows  that  the  largest  changes  in  velocity  occur  near  the 
rotor  inner  wall.  At  a position  of  Y/R^  = 0.04  near  the  rotor  inner  wall, 
operating  10%  low  in  flow  coefficient  (Basic  Flow  No.  2)  gives  an  increase 
of  12%  in  velocity  from  Basic  Flow  No.  1,  adding  a screen  (Basic  Flow  No. 
3)  gives  a decrease  of  5%  in  velocity  from  Basic  Flow  No.  1,  and  adding 
struts  (Basic  Flow  No.  4)  gives  an  increase  of  10%  in  velocity  from  Basic 
Flow  No.  1.  At  a position  of  Y/R^  = 0.40  near  the  rotor  mid-radius,  Basic 
Flow  No.  2 shows  a 4%  increase,  Basic  Flow  No.  3 shows  a 1%  decrease,  and 
Basic  Flow  No.  4 shows  a 2%  decrease  in  velocity  from  that  for  Basic  Flow 
No.  1. 

A comparison  between  the  inlet  velocity  data  and  the  critical 
cavitation  data  given  in  Figures  42  through  46  shows  large  changes  in 
cavitation  number  for  small  inlet  velocity  changes.  After  taking  the 
average  of  the  cavitation  data  at  each  flow  condition,  an  average  cavita- 
tion number  change  per  configuration  can  be  made.  This  change  gives  the 
mean  effect  of  variation  in  inlet  flow.  The  mean  effect  of  operating  10% 
low  in  flow  coefficient  is  to  increase  the  cavitation  number  by  46%,  the 
mean  effect  of  adding  a screen  is  to  decrease  the  cavitation  index  by  10%, 
and  the  mean  effect  of  adding  struts  is  to  increase  the  cavitation  number 
by  61%.  The  changes  in  inlet  velocity  near  the  inner  wall  are  less  than 
12%  as  discussed  in  the  previous  paragraph. 

The  measurements  made  near  the  end  of  the  rotor  cap  show  that  the 
structure  of  the  vortex  is  formed  along  the  inner  wall.  No  single 
experimental  curve  clearly  defines  how  this  vortex  developes.  However, 
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comparisons  of  experimental  data  do  give  insight  into  the  formation  of 
the  vortex. 

If  one  assumes  that  all  of  the  vorticity  associated  with  the  vortex 

is  in  the  vortex  core,  the  C associated  with  this  vortex  is  propor- 

^min 

tional  to  the  square  of  the  maximum  tangential  velocity  located  at  the 
core  radius  (r^)  of  the  vortex.  A comparison  of  Basic  Flow  No.  3 shown 
in  Figure  30  to  Basic  Flow  No.  1 shown  in  Figure  28  gives  a decrease  of 
12%  in  absolute  pressure  coefficient.  This  compares  favorably  to  a 
decrease  of  10%  in  cavitation  number.  However,  this  comparison  is  the 
only  one  which  correlates  with  cavitation  number.  As  an  example,  a 
comparison  between  Basic  Flow  No.  4 shown  in  Figure  31  and  Basic  Flow  No. 
1 shown  in  Figure  28  shows  a decrease  in  velocity  at  the  core  radii  but 
the  corresponding  cavitation  data  give  a 61%  increase.  In  addition,  the 
comparisons  between  the  measured  pressure  coefficients  at  the  end  of  the 
cap  and  the  average  cavitation  data  were  poor.  This  is  shown  in  Figures 
47  and  48.  The  lack  of  direct  correlation  of  cavitation  data  with 
measured  velocity  profiles  and  pressure  coefficients  implies  that  the 
minimum  pressure  of  the  vortex  does  not  occur  on  the  surface  of  the  rotor 
cap. 

A correlation  between  the  change  in  inlet  velocity  profile  and 
cavitation  data  and  a lack  of  correlation  between  the  measured  vortex 
tangential  velocity  and  cavitation  data  does  indicate  that  secondary 
flows  are  important.  Secondary  flows  generate  additional  streamwise 
vorticity  near  the  inner  wall  when  a boundary  layer  flow  is  turned  by  a 
rotor.  This  vorticity  will  induce  an  additional  component  of  Langentia] 
velocity  downstream  of  the  rotor  plane.  Whether  or  not  the  influence  of 
this  additional  streamwise  vorticity  on  the  tangential  velocity  is 
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completed  at  the  end  of  the  rotor  cap  is  not  important.  The  important 
fact  is  that  this  vorticity  will  change  the  net  circulation  near  the 
inner  wall. 

If  this  secondary  vorticity  is  an  important  part  of  the  net  circula- 
tion associated  with  the  vortex,  then  the  vortex  core  measured  at  the  end 
of  the  rotor  cap  should  be  proportional  to  the  distance  in  the  normal 
direction  from  the  surface  over  which  the  inlet  velocity  profile  is 
important.  In  addition,  the  vortex  core  should  be  inversely  proportional 
to  the  minimum  pressure  coefficient.  The  average  measured  vortex  size 
can  be  obtained  from  Figures  28  through  34.  Operating  10%  low  in  flow 
coefficient  (Basic  Flow  No.  2)  gives  a decrease  of  28%  in  core  radius 
from  Basic  Flow  No.  1,  adding  a screen  (Basic  Flow  No.  3)  gives  an 
increase  of  5%  in  core  radius  from  Basic  Flow  No.  1,  and  adding  struts 
(Basic  Flow  No.  4)  gives  a decrease  of  42%  in  core  radius  from  Basic  Flow 
No.  1.  These  results  correlate  the  trends  in  cavitation  data. 

There  are  two  possible  flow  parameters  which  control  the  core  radius 
of  the  vortex.  The  first  of  these  is  a boundary  layer  parameter,  such  as 
momentum  thickness.  The  other  one  ir  ^he  amount  of  passage  vorticity 
which  exists  near  the  inner  wail  of  the  rotor. 

The  boundary  layer  results  for  the  rotor  inlet  velocity  data  are 
given  in  Table  2.  These  calculations  show  that  operating  10%  low  in  flow 
coefficient  gives  a decrease  of  8%  in  momentum  thickness  from  Basic  Flow 
No.  1,  adding  a screen  gives  an  increase  of  15%  momentum  thickness,  and 
adding  struts  gives  an  increase  of  8%  in  momentum  thickness.  These 


results  do  not  correlate  the  trends  in  cavitation  data. 


4.  THEORETICAL  ANALYSIS  OF  FLOW  FIELD 


4 . 1 Introduction 

The  experimental  investigation  shows  that  the  structure  of  the 
vortex  is  formed  along  the  inner  wall.  Two  major  conclusions  can  be  made 
about  the  vortex.  Firstly,  the  critical  cavitation  number  of  the  vortex 
depends  on  the  inlet  velocity  gradient  near  the  inner  wall.  Secondly, 
the  vortex  core  radius  does  not  directly  correlate  with  any  boundary  layer 
parameter  of  the  inlet  velocity  profile. 

In  addition,  the  lack  of  correlation  between  the  measured  minimum 
pressure  coefficients  on  the  surface  of  the  rotor  cap  and  the  cavitation 
data  indicates  that  either  the  true  minimum  pressure  associated  with  the 
vortex  is  not  on  the  rotor  cap  or  vortex  breakdown  occurred  before  the 
end  of  the  rotor  cap.  In  either  case,  it  appears  that  the  flow  at  the 
rotor  exit  is  a better  indicator  of  the  vortex  behavior  than  at  the  end 
of  the  rotor  cap. 

The  calculation  of  the  flow  field  near  the  inner  wall  of  the  rotor 
is  a difficult  problem.  Assuming  that  the  flow  does  not  separate  along 
the  inner  wall,  the  flow  field  near  the  inner  wall  depends  upon  the 
potential  flow  effects  and  viscous  flow  effects  with  the  addition  of  the 
blade  to  blade  secondary  vorticity  effects.  For  simplicity  in  the  calcu- 
lation of  the  flow  field,  the  potential  flow  effects  are  the  axisymmetr ic 
flow  analysis  obtained  from  the  streamline  curvature  equations.  This  is 
referred  to  as  the  primary  flow.  The  blade  to  blade  secondary  vorticity 
analysis  is  referred  to  as  the  secondary  flow.  Although  viscous  effect  ; 
generate  the  inlet  velocity  profile  to  the  rotor,  no  viscous  effects  are 
included  in  the  calculation  of  the  primary  or  secondary  flow  through  the 
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rotor.  The  final  flow  field  analysis  through  the  rotor  requires  a 
coupling  of  the  primary  and  secondary  flows. 

In  order  to  calculate  the  secondary  flow,  a reasonable  estimate  of 
the  streamlines  near  the  inner  wall  is  necessary.  An  approximate  method 
lias  been  developed  to  calculate  these  mean  streamlines.  This  flow  field 
is  assumed  to  be  axisymmetric , inviscid,  and  incompressible.  The  boundary 
conditions  for  the  primary  flow  calculation  are:  (1)  the  flow  outlet 
angles,  (2)  the  velocity  profile  far  upstream  of  the  rotor  plane,  and  (.1) 
the  bounding  streamlines  of  the  flow.  The  flow  outlet  angles  are  obtained 

i, 

from  the  blade  angles  of  the  rotor  and  a deviation  angle.  A schematic'of 
the  boundary  conditions  is  shown  in  Figure  49. 

This  last  boundary  condition  is  the  most  difficult  to  define.  The 
open  rotor  is  unique  because  of  the  uncertainty  of  the  bounding  stream 
surface  of  the  flow  field.  The  flow  field  lacks  a well-defined  boundary, 
such  as  the  walls  of  a turbine  or  compressor.  For  this  reason,  the  flow 
field  for  the  open  rotor  is  considered  to  be  a region  constrained  within 
a potential  flow  stream  tube.  The  bounding  stream  surface  of  this  stream 
tube  is  actually  located  infinitely  far  from  the  rotor;  however,  it  is 
assumed  that  there  is  a finite  distance  beyond  which  the  flow  field  is 
relatively  unaffected  by  the  presence  of  the  rotor.  This  distance  is 
assumed  to  be  at  least  eight  times  the  rotor  radius.  This  stream  surface 
located  by  this  distance  will  serve  as  a flow  boundary.  In  this  manner, 
the  bounding  stream  tube  becomes  analogous  to  the  compressor  and  turbine 
walls . 

It  is  important  to  realize  that  these  boundary  conditions  imply  a 
solution  of  the  direct  turbomachinery  problem  and  not  the  indirect  or 
design  turbomachinery  problem.  For  the  direct  problem,  a rotor 
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configuration  is  specified  and  the  performance  for  any  given  initial 
condition  is  to  be  calculated.  Thus,  the  flow  outlet  angles  are  specified 
instead  of  the  tangential  velocity  (V  ) . 

The  equations  which  are  commonly  applied  to  calculate  the  flow  field 
in  an  axisymmetric  turbomachine  are  solved  by  use  of  the  so-called  stream- 
line curvature  method.  This  method  has  been  used  for  some  time  for  the 
design  and  analysis  of  various  types  of  turbomachines.  As  described  by 
Novak  l_  34_7,  Frost  j_  35_/,  and  Davis  j_  36_7,  this  system  of  equations  has 
been  successfully  applied  for  the  indirect  analysis  of  compressors  and 

turbines.  Also,  McBride  f_  37 / has  applied  this  streamline  curvature 

method  to  the  indirect  problem  of  an  open  hydrodynamic  rotor. 

In  this  chapter,  a method  to  calculate  the  final  mean  streamlines  for 
the  direct  problem  is  discussed.  The  method  uses  the  streamline  curvature 
equations  to  establish  the  flow  field  from  the  boundary  conditions.  The 
equations  used  for  secondary  flow  calculations  are  discussed. 

4 . 2 Primary  Flow  Field 

4.2.1  Discussion  of  Calculation  Procedure.  A schematic  of  the 
calculation  procedure  for  the  final  flow  field  is  given  in  Figure  50. 

This  outlines  the  iterative  procedure  for  the  calculation  of  the  final 
flow  field  and  indicates  where  secondary  flow  calculations  are  employed. 

The  initial  conditions  (Step  1)  to  the  solution  of  the  direct  problem 
are  (1)  bounding  stream  tube  and  (2)  velocity  profile  in  rotor  plane  with- 
out rotor.  With  this  information,  the  initial  streamlines  without  rotor 
can  be  calculated  (Step  2).  The  result  of  this  calculation  is  the 
boundary  condition  of  an  initial  velocity-energy  profile  at  a station  far 
upstream  of  the  rotor  plane. 
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Figure  50.  Outline  of  Calculation  Procedure 
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Knowing  the  blade  angles,  the  first  estimate  of  the  flow  outlet 

angles  (Step  3)  can  be  calculated.  These  flow  outlet  angles  depend  on 

the  blade  angles  (3  ) and  deviation  angle  (6).  Hie  deviation  angle 
2 

correlation  developed  by  Howell  / 38/  is  initially  applied.  This  rela- 
tionship considers  only  thin  blade  sections  and  assumes  that  each,  blade 
section  operates  near  design  incidence.  As  shown  in  Figure  49,  all  of 
the  boundary  conditions  are  now  known  and  the  flow  field  can  be  solved 
with  the  rotor  included  (Step  4) . 

Once  a convergent  solution  is  obtained  for  the  flow  field  using 
Howell’s  deviation  rule  (Step  4),  the  axial  velocity  distributions  are 
known  whereby  the  inlet  angles  can  be  estimated  in  addition  to  the  accel- 
eration through  the  rotor.  Now  a second  estimate  of  the  flow  outlet 
angles  (Step  5)  can  be  made.  For  this  deviation  angle,  the  effects  of 

■k 

acceleration  (AS'),  blade  camber  (S  ),  and  blade  thickness  (AS  ) are 

o 

calculated  separately.  For  the  calculation  of  the  deviation  term  due  to 

axial  acceleration  through  the  rotor,  an  equation  developed  by 

Lakshminarayana  [_  39 / is  applied.  For  the  calculation  of  deviation  terms 

due  to  camber  and  thickness  effects,  the  data  collected  by  Lieblein  j_  40 / 

are  used.  The  result  is  an  improved  outlet  flow  angle  profile  which  can 

be  used  to  again  calculate  the  flow  field  (Step  6). 

The  convergent  solution  of  the  flow  field  (Step  6)  is  then  used  to 

solve  the  secondary  vorticity  equations  (Step  7)  and  to  determine  a 

deviation  term  (A5  ) which  is  due  to  blade  to  blade  secondary  flows.  The 
s 

details  of  the  secondary  flow  calculations  are  described  in  Section  4.3. 

An  improved  outlet  flow  angle  profile  (Step  8)  is  obtained  by  adding  this 
secondary  flow  term  to  the  deviation  terms  previously  calculated  to  obtain 


| 
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3 * = 6 - A5  * + A6*  + 6 + A6  , (21) 

2 2 os 

where  3 is  the  outlet  flow  angle  and  S is  the  blade  outlet  angle. 

2 2 

Tilts  outlet  flow  angle  profile  is  then  used  as  a boundary  condition  In 
the  calculation  of  the  final  flow  field  (Step  9) . 

Finally,  all  of  the  deviation  angle  calculations  are  checked  based 
on  the  flow  field  calculated  in  Step  9.  If  the  angles  did  not  change 
significantly , then  the  results  obtained  in  Step  9 are  the  final  flow 
field.  If  the  angles  were  different,  then  these  new  angles  are  used  to 
again  calculate  the  final  flow  field  (Step  9) . 

4.2.2  Streamline  Curvature  Method.  The  major  equations  used  in  the 
streamline  curvature  method  of  anal'  si  are  derived  from  the  principles 
of  conservation  of  mass,  momentum,  and  total  energy.  In  this  analysis, 
the  fluid  is  assumed  to  be  incompressible,  inviscid,  and  steady.  Because 
the  flow  field  is  axisymmetric , it  is  necessary  to  use  only  two  components 
of  velocity.  One  component  of  velocity  (the  meridional  velocity  (V  ) ) is 
tangent  to  the  streamline  and  is  projected  onto  the  meridional  plane. 

This  component  is  related  to  the  axial  velocity  by  the  cosine  of  the 
streamline  angle  (4>)  . The  second  component  of  velocity  is  called  the 
tangential  velocity  (V  ) and  is  in  the  circumferential  direction  which  is 
perpendicular  to  the  meridional  plane. 

The  resultant  equations  allow  for  streamline  curvature  and  for 
vorticity  in  the  flow.  However,  It  is  important  to  realize  that  the 
solution  to  the  flow  field  does  not  contain  all  of  the  vorticity.  In 
particular,  only  the  circumferential  vorticity  is  totally  Included.  The 
other  components  of  vorticity  contain  derivatives  with  respect  to  the 
circumferential  direction  which  are  assumed  to  be  zero  with  an 
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axisymmetric  flow  analysis.  As  discussed  by  Hawthorne  and  Novak  l_  19 
the  neglected  vorticity  terms  can  be  related  to  the  secondary  flows  that 
occur  in  the  blade  passage. 

The  set  of  equations  used  for  this  streamline  curvature  method  are 
as  follows: 

r-momentum  equation  — 


1 3JP 
p 3r 


cosft 

R 

m 


sln<j>  m 

V 3m 
m 


9-momentum  equation  — 


3(ry 

39 


= 0 


(22) 


(23) 


energy  equation  — 

3P 

- — = 0 {except  across  rotor  plane)  , (24) 


and  continuity  equation  — 


r 

(2 


2ir 


V cosi))  r dr  = constant 
m 


(25) 


In  these  equations,  R is  the  radius  of  curvature  of  the  streamlines  in 

m 

the  meridional  plane  and  <J>  is  the  streamline  angle  with  respect  to  the 
axial  direction. 

In  Equation  (22) , there  are  three  contributing  factors  in  the 
development  of  the  radial  pressure  gradient.  The  first  term  is  directly 
related  to  the  centrifugal  force  which  the  fluid  experiences  upon  passing 
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through  the  rotor.  The  second  term  is  due  to  a moving  fluid  particle 
being  subjected  to  the  streamline  curvature  in  the  meridional  plane.  The 
magnitude  of  this  term  depends  directly  on  the  radius  of  curvature  of  the 
streamline.  The  third  term  is  due  to  the  convective  acceleration  as  the 
flow  area  either  converges  or  diverges.  The  combination  of  these  three 
terms  defines  the  radial  pressure  gradient. 

The  computational  procedure  using  the  four  equations  is  an  iterative 
one.  Initially,  the  velocity  distribution  at  the  reference  station  far 
upstream  of  rotor  is  transferred  to  the  downstream  flow  by  using  the 
continuity  equation.  At  the  rotor  exit  plane,  the  tangential  velocity 
profile  is  related  to  the  axial  velocity  by 

V0  = ilr  - Vx  tanS^*  , (26) 

2 2 

* 

where  2 is  the  flow  outlet  angle.  Now,  the  principle  of  constant 
2 

angular  momentum  allows  for  the  transfer  of  tangential  velocity  downst re... : 
of  the  rotor.  Finally,  the  static  pressure  can  be  calculated  along  the 
streamlines  using  the  radial  pressure  equation  and  a new  meridional 
velocity  profile  is  calculated  by  using  this  pressure.  This  velocity- 
pressure  coupling  is  repeated  until  the  streamlines  are  located  such  that 
the  flow  simultaneously  satisfies  continuity,  conservation  of  momentum, 
and  conservation  of  total  energy. 

4.2.3  Estimation  of  Rotor  Outlet  Angles.  The  accuracy  of  the 

calculated  velocity  profiles  to  experimental  data  depends  primarily  upon 

* 

the  prediction  of  the  flow  outlet  angles  (B  ) used  in  Equation  (26). 

2 

* 

The  flow  outlet  angle  (8  ) is  shown  in  Figure  51.  This  angle  is  a 

2 

combination  of  blade  geometry  and  a deviation  angle.  The  blade  angle 


Figure  51.  Definition  of  Blade  Angles 
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(3  ) is  defined  as  the  angle  between  the  tangent  to  the  blade  section 
camber  line  at  the  trailing  edge  and  the  exit  axial  direction. 

The  correct  determination  of  the  deviation  angle  presents  a problem. 
In  particular,  for  the  case  of  a rotor  having  a low  hub-tip  ratio  and 
twisted  blades,  the  deviation  angle  is  more  difficult  to  correlate  with 
blade  geometry  even  though  the  flow  is  turned  through  a small  angle. 
However,  the  flow  turning  is  still  an  expression  of  the  guidance  of  the 
passage  formed  by  adjacent  blade  sections,  and  it  is  expected  that  the 
geometry,  such  as  blade  camber  (0  ) , thickness  (t),  solidity  (c/s),  and 
stagger  angle  (1),  be  the  principle  influencing  factors.  Thus,  deviation 
angles  can  be  estimated  for  each  rotor  radius  based  on  the  geometry  of 
the  blade.  It  can  be  expected  that  the  largest  errors  will  occur  at  the 
blade  tip  where  three-dimensional  factors  dominate  and  at  the  inner  wall 
where  viscous  effects  dominate. 

From  extensive  cascade  potential-flow  theories,  it  was  found  that 

t 

the  deviation  angle  increases  with  blade  camber  and  chord  angle  and 
decreases  with  solidity.  Using  this  result,  many  experimentors  have 
correlated  data  to  determine  a deviation  rule.  In  particular,  Howell 

j_  38 / has  developed  a correlation  on  the  basis  of  nominal  operating 

conditions.  This  condition  pertains  to  a deflection  of  the  flow  which 
is  80%  of  its  maximum  stalling  deflection.  The  deviation  angle  is  defined 
for  each  blade  section  as 

0.23(9  ) (s/c) 1 /2 (2a/ c) 2 
6 = £ 

H 1.0  - (0.1/50X9  Xs/c)1/2 

c 


(27) 
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where  6 is  the  section  turning  angle,  s/c  is  the  space-to-chord  ratio, 
c 

and  a/c  is  the  distance  from  the  section  leading  edge  to  the  point  of 
maximum  camber  divided  by  the  section  chord  length. 

Equation  (27)  is  used  in  the  first  estimate  of  flow  outlet  angles 
(Step  4)  in  the  calculation  of  the  flow  field.  It  is  important  to  note 
that  this  deviation  correlation  assumes  thin  blades  and  that  each  blade 
section  operates  near  design  incidence.  Once  a converged  solution  is 
obtained  for  the  flow  field  using  Howell's  deviation  (Step  5),  the  inlet 
and  outlet  axial  velocity  distributions  are  estimated.  Now,  the  deviation 
angle  can  be  recalculated  based  on  the  effects  of  axial  acceleration 

•k 

(AS  ),  camber  (S^) , blade  thickness  (AS  ),  and  secondary  flows  (A5g) . 
Therefore,  a second  order  deviation  angle  can  be  expressed  as 

5 = 5 + AS*  - AS'  + AS  . (28) 

o s 


2irk(c/s) tang 

1 
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+ tan0 


9 n 

cos  8 


(29) 


where  AVR  is  the  axial  velocity  ratio  (V  /V  ) , k is  the  cascade 

x x 
2 1 

influence  coefficient,  G/c  is  the  ratio  of  the  distance  from  the  chord 
line  to  the  maximum  camber  point  of  the  section  to  chord  length,  6 and 

•k  * 

6 are  the  inlet  and  outlet  flow  angles,  respectively,  and  a is  defined 

2 

as  the  difference  between  the  iniet  flow  angle  (3  ) and  the  section 
stagger  angle  (X) . This  relationship  is  valid  for  small  changes  in  AVR 
and  for  small  turning  angles.  Perhaps  most  important,  the  deviation  angle 
due  to  axial  acceleration  in  the  flow  always  reduces  the  total  deviation. 
The  deviation  due  to  thickness  and  camber  is  given  by  empirical  data 

which  have  been  collected  by  Lieblein  J_  40 / . The  thickness  deviation  is 

calculated  from 


A6*  = <Vsh  • (Vt  • <O10  • {30) 

where  (6°)  represents  the  basic  variation  for  the  10%  thick  65-series 

o 10 

thickness  distribution,  (K.)  represents  any  correction  necessary  for  a 

o sh 

blade  shape  with  a thickness  distribution  different  from  that  of  the 

65-series  blade,  and  (Kj)  represents  any  correction  necessary  for  maximum 

blade  thickness  other  that.  10%.  The  value  of  (K,)  can  range  from  1.1  to 

o sh 

0.7;  therefore,  an  average  value  of  0.9  is  used  for  the  rotor.  Also,  the 
other  coefficients  are  determined  from  graphs  of  experimental  data  in 
Reference  ]_  4 0_/  knowing  solidity  and  inlet  flow  angles. 
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The  effect  of  camber  on  deviation  is  calculated  by  using 

5 = m 9 , (31) 

o c 

where  m is  the  slope  of  the  deviation  angle  variation  with  camber  angle 
(9c) . Values  of  m can  be  determined  from  the  figures  given  in  Lieblein 

j_  40 / by  knowing  the  flow  inlet  angle  and  solidity.  This  slope  is  given 

for  either  circular-arc  or  parabolic-arc  mean  camber  lines. 

Equations  (29),  (30),  and  (31)  are  used  for  the  second  estimate  of 
rotor  outlet  angles  (Step  5) . In  order  to  use  these  deviation  relation- 
ships the  axial  velocity  through  the  rotor  had  to  be  estimated.  This 
velocity  estimation  was  accomplished  by  using  Howell's  correlation  in  the 
first  calculation  of  the  flow  field. 

The  last  deviation  term  to  be  considered  is  due  to  secondary  flows. 
The  effects  of  secondary  flows  are  most  dramatic  near  the  rotor  inner  wall 
where  the  inlet  velocity  gradients  are  largest.  Secondary  flows  produce 
a streamwise  component  of  vorticity  which  is  not  taken  into  consideration 
by  axisyrametric  analysis. 

The  solution  for  the  streamwise  secondary  vorticity  involves  the 
numerical  solution  of  two  equations  which  are  discussed  in  the  next 
section.  After  solving  for  the  streamwise  vorticity  at  the  rotor  exit 
plane,  a secondary  stream  function  is  defined  in  this  plane.  The  solution 
for  the  stream  function  gives  perturbation  velocities  which  are  used  to 
find  a deviation  angle  (A6g) . 

4 . 3 Secondary  Flow  Field 

4.3.1  Secondary  Vorticity.  In  order  to  accomplish  work  in  a rotor, 
angular  momentum  must  be  imparted  to  the  fluid.  The  lifting  action  cf  ' e 


A 
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rotor  blades  produces  tangential  pressure  gradients.  For  the  case  of  the 
rotor  operating  within  a boundary  layer,  the  pressure  gradients  are 
balanced  by  streamline  curvature  at  a sufficient  distance  from  the  inner 
wall.  Close  to  the  inner  wall,  the  velocities  are  much  less  than  at  the 
rotor  tip.  The  fluid  near  the  inner  wall  does  not  have  sufficient 
momentum  to  balance  the  pressure  gradients  imposed  by  the  flow.  The 
result  is  a tangential  component  containing  vorticitv  aligned  in  the 
streamwise  direction.  This  additional  streamwise  vorticity  causes  a 
deviation  in  the  outlet  angles. 

There  are  no  simple  techniques  that  accurately  predict  these  inner 
wall  velocity  gradient  effects  for  rotational  flows.  It  is  important  to 
recognize  that  a satisfactory  description  of  this  boundary  layer  flow 
cannot  come  from  refinements  of  the  two-dimensional  boundary  layer  theory, 
but  rather  from  three-dimensional  analysis.  The  reason  is  that  boundary 
layer  behavior  in  a rotor  exhibits  variations  in  its  tangential  direction. 

The  occurrence  of  these  secondary  flows  which  are  boundary  layer  flows 
having  a component  normal  to  the  mainstream  direction  arise  principally 
from  blade  to  blade  and  radial  pressure  gradients.  Similarity  between  the 
results  of  secondary  flow  analysis  and  the  cross-flow  in  the  outer  part  of 
a three-dimensional  boundary  layer  has  been  established  by  Horlock  J_  41_/. 

The  distributed  passage  secondary  vorticity  generated  by  turning  a 
boundary  layer  flow  is  one  of  three  possible  types  of  secondary  streamwise 
vorticity,  which  can  occur  near  the  inner  wall  of  a rotor.  As  discussed 
by  Hawthorne  / 42_/,  three  components  of  secondary  vorticity  are  identi- 


fied in  the  direction  of  flow  at  the  exit  of  a blade  row.  These  secondary 
sources  of  vorticity  are  usually  regarded  as  a perturbation  on  the  primary 
flow.  They  are  (1)  a distributed  passage  vorticity  in  the  blade  passage 
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which  may  result  in  the  formation  of  the  so-called  passage  vortex,  (2) 
the  trailing  shed  vorticity,  and  (3)  the  trailing  filament  vorticity. 

The  latter  two  types  of  vorticity  are  due  to  the  passage  vortex  and  occur 
in  the  vortex  sheet  leaving  the  blade  trailing  edge.  They  lead  to  the 
formation  of  another  vortex  which  is  opposite  in  rotation  to  that  of  the 
passage  vortex  as  shown  in  Figure  52.  The  trailing  filament  vorticity  is 
caused  by  the  stretching  of  the  vortex  filaments  as  they  move  over  the 
surface  of  the  blades.  The  last  component  of  vorticity  is  the  trailing 
shed  vorticity  which  is  caused  by  a change  in  circulation  along  the  span 
of  the  blades. 

Each  of  these  vorticity  components  could  be  attributed  to  the  exist- 
ence of  the  wall  boundary  layer  and  not  the  change  of  the  boundary  layer 
due  to  viscous  effects  as  flow  passes  through  the  rotor.  Therefore,  the 
primary  assumption  leading  to  the  existing  theoretical  descriptions  of 
secondary  flows  is  that  viscous  effects  produce  a boundary  layer  on  the 
inner  wall  upstream  of  the  rotor.  Whereas  within  the  rotor,  the  imposed 
pressure  gradients  play  the  major  role  and  viscosity  has  little  effect  in 
the  resulting  secondary  flows.  This  assumption  is  characteristic  of  what 
is  generally  termed  Inviscid  secondary  flow  analysis. 

As  shown  in  an  analysis  of  secondary  flows  by  Came  and  Marsh  t_  43/, 
the  total  strength  of  these  three  secondary  components  of  streamwise 
vorticity  is  zero  in  the  flow  downstream  of  a many  bladed  cascade. 
However,  this  vorticity  does  have  an  effect  on  the  flow  field.  The 
primary  effect  of  this  secondary  vorticity  is  to  cause  a deviation  in  the 
flow  outlet  angle.  This  is  due  to  the  induced  velocity  field  from  the 
passage  secondary  vorticity.  This  deviation  can  be  quite  large  when  the 
inlet  velocity  gradient  to  the  rotor  is  large  or  when  the  flow  is  turned 
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through  a large  angle.  The  solution  for  this  deviation  must  be  consistent 
with  the  trailing  vortex  sheet,  but  the  strength  of  this  sheet  does  not 
need  to  be  known. 

Theoretical  development  of  secondary  flow  equations  begin  with  the 
fluid  flow  equations  that  govern  the  flow  process.  The  equation  of  motion 
for  incompressible  flow  with  reference  to  axis  rotating  at  constant 
angular  velocity  (SI)  is  given  in  Greenspan  _/  4 4 / as 

TxW  + 20xW  = - V*  (")  - V'[w2/2  - 1/2 (f2xr) 

• (iTxrj]  - (u/p)V'xV'xW  (32) 

where  the  prime  denotes  differentiation  with  respect  to  the  rotating  frame 
(SI) , W is  the  relative  velocity,  and  C is  the  relative  verticitv  giver,  by 
\ = V'xW  = u - 2fi.  Taking  the  curl  of  Equation  (32)  gives  the  vorticity 
eauation 


(W-V')T  = (T*V')W  - 2V 'x(Hx"w)  + (p/o)V'21  . (33) 


Expressions  for  the  absolute  secondary  vorticity  (u)  defined  along  a 
rela  ive  streamline  valid  for  a rotor  are  obtained  from  Equation  (33)  by 
Lakshminarayana  and  Horlock  j_  16 / . The  resultant  equations  are 


2u  , 2n,  ,0)  , 

n b n 

* o 

WR  M2 


n ' b ' + 

W2 


(p/o)s' 


(34) 


and 


3 i-  Wub ' 

an7"  L“>n.wJ  = ~ 


Woo 


8 a 


b' 


*b' 


3s' 


+ (u/o)n'  • V’2C 


(35) 


102 


where  as  shown  in  Figure  53,  s',  n',  b'  represent  the  natural  coordinates 
for  the  relative  flow,  W is  the  relative  velocity,  oj^ , and  are 

absolute  vorticitv  resolved  along  the  relative  streamline  direction  (s') 


and  the  principal  normal  direction  (n1),  Ci  is  the  rotation  vector,  and  R 
is  the  radius  of  curvature  of  the  relative  streamline.  In  this  relative 
coordinate  system  indicated  by  the  primes,  the  streamwise  direction  (s') 
is  defined  along  the  flow  direction,  the  principal  normal  direction  (n ' ) 
is  defined  toward  the  center  of  curvature  of  the  streamlines,  and  the 
bi-normal  direction  (b’)  is  defined  as  s’xn'  = b'.  Additional  definitions 
of  the  flow  in  the  relative  coordinate  system  are  given  in  Reference 


L i6_/- 

The  means  by  which  the  streamwise  component  of  vorticity  is  produced 
in  a relative  flow  are  similar  to  those  in  a stationary  system.  However, 
it  is  important  to  note  that  additional  secondary  vorticity  is  generated 
when  flxW  has  a component  in  the  relative  streamwise  direction.  Rotation 
has  no  effect  when  the  absolute  vorticity  vector  lies  in  the  s'-n'  plane 
and  the  rotation  (f.)  has  no  component  in  the  bi-normal  direction  (b1). 

The  equations  for  secondary  vorticity  created  in  the  passage,  namely. 
Equations  (34)  and  (35),  are  extremely  difficult  to  apply  because  they  are 
nonlinear  partial  differential  equations.  Therefore,  some  simplifying 
assumptions  and  specif ications  of  the  flow  field  must  be  made  for  this 
application.  The  basic  flow  parameters  for  this  vortex  problem  are  a 
thick  boundary  layer  entering  the  rotor  and  the  rotor  operating  at  a flow 
coefficient  less  than  0.5.  It  Is  also  important  to  note  that  the  primary 
flow  through  the  rotor  is  rotational. 

In  order  to  evaluate  the  relative  importance  of  the  terms  that 
describe  the  passage  secondary  vorticity,  the  vorticity  equations  were 
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first  nondimensionalized.  Then  an  order  to  magnitude  evaluation  was 
applied.  The  details  of  this  process  are  given  in  Reference  / 45 _/. 
The  resulting  normal  relative  vorticity  equation  is 


3 

3s  f 


(<un.W) 


Wm  , 
n 

V 


(36) 


where  a , is  the  distance  between  neighboring  streamlines  in  b'  direc- 
b 

tion.  Two  assumptions  applied  in  developing  Equation  (36)  are  (1)  the 
ratio  of  boundary  layer  thickness  to  rotor  tip  radius  is  of  order  one  and 
(2)  the  blade  chord  (c)  is  an  order  of  magnitude  less  than  the  camber 
radius . 

Equation  (36)  can  be  integrated  along  a relative  streamline  by 
combining  the  two  terms  into  one  differential.  The  result  is 


w 


1 


W a,  , 

1 


w v 


(37) 


where  the  subscript  1 refers  to  the  rotor  inlet  and  unnumbered  subscripts 

refer  to  any  position  along  a streamline  in  the  rotor.  This  simplified 

equation  shows  a dependence  of  the  relative  component  of  the  normal 

vorticity  on  the  variation  of  relative  velocity  (W)  along  a relative 

streamline,  on  the  relative  component  of  the  absolute  inlet  normal 

vorticity  (oi^,),  and  on  the  flow  convergence-divergence  in  the  bi-normal 
1 

direction  (a,  ,). 

b 

Equation  (37)  can  be  evaluated  along  a relative  streamline  if  the 
streamline  spacing  (a  , ) in  the  bi-normal  direction  and  the  relative 

D 

velocity  along  a relative  streamline  are  known.  These  quantities  can  be 
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obtained  from  the  primary  flow  analysis  using  the  streamline  curvature 
equations . 

For  some  rotors  having  a high  hub-to-tip  ratio,  the  streamline 
spacing  in  the  radial  direction  (a^,)  remains  constant,  so  that  Equation 
(37)  becomes 


1 


(o  , = u>  , • rr 

n n W 
1 


(38) 


This  result  relates  the  change  in  the  relative  component  of  the  absolute 
normal  vorticity  (u)^,)  t0  changes  in  the  relative  velocity  (W) . Also, 
assuming  a weak  shear  flow  so  that  the  axial  velocity  remains  constant 
gives 


u , = 0)  , 
n n 


cosB 


cosB 


(39) 


These  last  two  relationships  are  commonly  used  for  secondary  flow  cascade 
theory. 

The  resulting  streamwise  vorticity  equation  is 


3s' 


“s' 


2(o  , 20  ,io  , 20  ,io  , 

> = n'  + b'  n _ n'  b 


(A0) 


WR 


Wz 


The  second  term  on  the  right-hand  side  of  Equation  (A0)  is  due  to  curva- 
ture induced  secondary  vorticity.  The  third  term  is  due  to  rotation 
induced  secondary  vorticity.  The  last  term  is  important  because  0 , is 
large  because  the  flow  coefficient  is  less  than  one.  Integrating 
Equation  (A0)  in  the  relative  streamwise  direction  yields 
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u>  , = W 
S2  2 


2ui  , 


WR 


ds'  + W 


2W 


ds' 
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n b , i 

* 

2 

ds  + to  . ■ 

,.o  S 

w 

4 

w-  , 

1 

1 1 

It  is  important  to  note  that  the  boundary  layer  thickness  Is  the 

same  order  as  the  rotor  tip  radius  and  much  larger  than  the  blade  chord. 

The  velocity  gradients  for  these  conditions  are  not  as  large  as  would  he 

the  case  if  the  boundary  layer  thickness  would  be  much  less  than  the 

rotor  tip  radius.  Therefore,  stream  surfaces  in  the  blade  passage  will 

not  be  significantly  distorted  so  that  Bernoulli's  surface  rotation  (see 

Reference  ]_  46 _/)  is  not  important. 

If  the  relative  streamlines  and  the  absolute  vortex  lines  lie  on 

cylindrical  surfaces  (to,  , = 0)  and  the  rotation  vector  (Q)  is  parallel  to 

b 

the  axis  of  the  cylindrical  surfaces,  then  Equation  (41)  would  become 


s s 

2 1 

W 

2 1 
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J WR* 
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ds' 


(42) 


This  relationship  is  commonly  used  in  most  cascade  flow  calculations  and 

is  the  equivalent  of  an  equation  developed  by  Squire  and  Winter  j_  4 7 / 

but  with  a relative  reference  frame. 

More  comparisons  between  the  equations  for  the  secondary  vorticity  in 
a rotating  reference  frame  and  cascade  secondary  vorticity  are  discussed 
in  Reference  / 45  /. 


The  solution  for  the  absolute  streamwise  vorticity  defined  along  a 
relative  streamline  can  be  calculated  using  Equations  (37)  and  (41) . 
However,  the  primary  flow  field  must  be  initially  calculated. 

From  the  solution  of  the  primary  flow  field,  the  relative  velocity 

(W) , the  streamline  spacing  in  the  b'-direction  (a^,),  the  rotation  (..) , 

* 

the  flow  turning  angle  (3  -3  ),  and  the  inlet  absolute  normal  velocity 

1 2 

(u>  ) are  known  for  a mean  relative  streamline, 

n 

1 

As  shown  in  Figure  54,  these  quantities  can  be  divided  into  their 

components  and  used  in  the  secondary  flow  equations.  A mean  radius  of 
* 

curvature  (R  ) for  the  blade  sections  defined  by  a mean  relative  stream- 
line is  calculated  from 


(6  -B  ) 

1 2 

for  small  turning  angles  (B  -g  ).  Also,  the  component  of  absolute 

1 2 

vorticity  defined  in  the  relative  bi-normal  direction  (w,  , ) is  assumed  to 

b 

be  proportional  to 


u 


b’ 


(44) 


for  the  rotor  considered. 

4.3.2  Secondary  Velocities.  Secondary  motions  are  induced  by  the 
distributed  secondary  vorticity  in  the  flow  between  the  blades.  With 
large  blade  turning  angles,  the  secondary  motions  in  the  blade  passage 
are  masked  by  large  spanwise  flows  associated  with  changes  in  axial 
velocity  profile  referred  to  as  displacement  effects.  Hawthorne  and 
Armstrong  J_  48 / have  developed  an  analysis  for  predicting  these 


T 
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displacement  effects  which  reduce  the  secondary  motions.  With  small 
blade  turning  angles,  the  secondary  motion  largely  determines  the  varia- 
tion in  outlet  angle  downstream.  For  most  rotor  design,  the  turning  near 
the  inner  wall  is  small  where  secondary  motions  are  the  largest. 

Formulas  for  calculating  the  two-dimensional  solution  for  the 
secondary  motion  in  blade  passages  have  been  given  by  Hawthorne  j_  49 _/, 
Hawthorne  and  Novak.  ]_  19 _/,  and  Dixon  l_  50 _/.  The  secondary  motion  in 
terms  of  velocities  (u,  v,  and  w)  is  assumed  to  be  induced  by  the  stream- 

wise  secondary  vorticity  (w  ,-ui  ,)  at  the  trailing  edge  plane  of  the 

s s 
2 1 

rotor.  This  streamwise  secondary  vorticity  was  calculated  along  a mean 
streamline  and  is  a function  of  radius  only.  It  is  important  to  note 
that  the  normal  component  of  vorticity  (w  , ) is  accounted  for  in  the 
nxi symmetric  flow  analysis. 

The  effect  of  the  streamwise  component  of  secondary  vorticity 

(u  , —co  ,)  within  the  blade  passage  is  similar  to  that  obtained  in  the 

2 1 _ _ _ _ 
flow  through  a duct  bend  (Hawthorne  j_  51  _/,  Eichenberger  j_  5 2 / , Detra 

]_  5 3 / ) ; however,  there  is  the  difficulty  of  devising  a reasonable 

approximate  method  of  satisfying  the  Kutta-Joukowski  condition  at  the 

trailing  edge. 

The  method  developed  by  Hawthorne  and  Novak  [_  19/  assumes  that  the 
flow  is  contained  in  a duct  defined  by  the  blades  and  streamlines  of  the 
primary  flow  leaving  the  trailing  edge  of  each  blade.  In  this  trailing 
edge  plane,  the  secondary  streamwise  vorticity  is  divided  into  a 

■h 

tangential  and  axial  component.  The  former  term,  namely,  (u  , -w  ,)sin£ 

2 1 2 

causes  a radial  gradient  of  axial  velocity  and  the  latter  leads  to  an 
equation  for  a stream  function  describing  the  radial  and  tangential 
velocities  in  the  trailing  edge  plane  (r,0). 


no 


Referring  to  Figure  55,  the  velocity  perturbations  induced  by 

(u  ,-u^,)  have  components  u,  v,  w in  the  r,  0,  x directions,  respectively. 
2 1 

The  secondary  streamwise  vorticity  as  calculated  from  Equations  (37)  and 
(41)  and  the  relative  velocity  (W)  are  only  a function  of  radius.  The 
components  of  the  streamwise  vorticity  are 


& riw 

(w  ,-u)  ,)sin8  = w = - -r-  (45) 

s s o 9 dr 

2 1 2 


and 


(u)  ,-u  ,)cos8  = w 
s s 2 x 

2 1 2 


i 3 (tv)  3 (u) 

r 9r  39 


(46) 


where  w is  only  a function  of  radius  and  u is  not  a function  of  x. 
Equation  (45)  can  be  integrated  to 


w 


(oj  ,-w  ,)sin8 

S S o 

2 1 


dr 


(47) 


where  the  constant  of  integration  is  obtained  from  continuity.  The  net 
mass  flow  through  the  exit  plane  due  to  the  secondary  velocity  (w)  must 
be  zero.  Therefore,  the  boundary  condition  for  Equation  (47)  is 


% 


2tt 


wr  dr  = 0 


(48) 
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Defining  a secondary  stream  function  ('F)  as 


1 3>F 

U=739  3nd 


jTF 

3r 


(40) 


the  continuity  equation  and  the  x-direction  vorticity  equation  can  be 
combined  into  a stream  function  equation.  The  result  is 


32V  1 9T  1_ 

3r2  r 3r  r2  392 


- (u  ,-u  t)cosB 
s s 2 

2 1 


(50) 


Equation  (50)  describes  the  radial  and  tangential  velocities  in  the 
trailing  edge  plane  (r,9).  The  boundary  conditions  are 

T = 0 at  r = Ri>RR  (51) 

and  at  the  blade  trailing  edge, 

* * 
w sinfi  = v cosB 
2 2 


or 


3Y  + a * 

— = - w tanB 
3r  2 


f(r) 


(52) 


at 


9 = 0 , 2tt/N  , WN 

Hawthorne  and  Novak  f_  19 / redefine  a new  stream  function  which 

includes  the  boundary  condition  of  Equation  (52).  The  result  is  that  the 
stream  function  will  now  be  zero  everywhere  on  the  boundary.  The  stream 
function  is  defined  as 


113 


k. 


Y( r,9)  = Y (r ,9)  + V (r) 


(53) 


where  t'  = 0 at  r = R^,R^  and  at  9 = 0,  2’r/N,  4tt/N, 
Eq>  .tion  (50)  becomes 


a2v*  , l 3**  ^ l a2**  w d , n *v 

+ - ^ — + = - — (r  tanB  ) 

r dr  ? r dr 
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dw 
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- (u  ,-(j  ,)cos8  . (54) 

S S o 

2 1 

Finally,  Equation  (45)  can  be  used  in  Equation  (54)  giving 

V2V  = — 3—  (r  tanB  ) - (oj  ,-w  ,)secS  . (55) 

r dr  o s s 2 

2 1 

The  right-hand  side  of  Equation  (55)  is  a function  of  r only.  Equation 

(55)  is  identical  to  that  obtained  by  Hawthorne  and  Novak  j_  19 / except 

the  equation  is  for  a relative  reference  frame. 

In  a discussion  of  a paper,  Hawthorne  j_  54  / gives  the  transforma- 
tion of  Stoke’s  equation  from  r,9  coordinates  to  z,y  coordinates.  The 
transformation  is 


r/R.  = e 
1 


(56a) 


and 


(56b) 


Substituting  Equation  (5b)  into  the  stream  function  equation  Rives 


+ Alt.  = (rJ-r)-'  F(r)  = (i(z) 

o o N 

<>y 


with  the  boundary  conditions  of  , = 0 at  y = 0 and  y = 1 and  at  z - 0 and 
z = it.  The  quantity  S.  is  defined  as 


N Ri 

1 * %>  ‘"t’ 
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The  solution  of  Equation  (57)  must  be  periodic  so  that 


G(z)  = G(z)  I ~r  sin  rnry 
1,3,5  n 


and  is  in  the  form 


T = Z V (z)  sin  m;y 

1,3,5  n 


A general  solution  of  Equation  (57)  in  the  form  of  Equation  (59)  has 
been  obtained  by  Hawthorne  MW  using  the  method  of  variation  of 
parameters.  The  solution  for  T (z) 


T (z)  = 
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and  for  V'(z)  is 
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A mean  deviation  angle  (A6^)  due  to  secondary  flows  can  b 
for  any  radius.  The  mean  deviation  angle  is  defined  as 


where 
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— * * 
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Using  Equations  (A9) , (52),  and  (53),  Equation  (62)  becomes 
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In  terms  of  the  transformed  z,v  coordinates,  Equation  (64)  can 


(61) 


>•  calcul  a ted 


(61) 


(63) 


(64) 


be  written 


as 


VORTEX  PARAMETERS 


5.  CALCULATION  OF  FLOW  FIELD  AND 

5 . 1 Calculat i a n of  Final  low  Fi eld 

5.1.1  Flow  Field  for  Basic  Flow  No.  1.  The  case  of  the  rotor 
operating  near  design  (Basic  Flow  No.  1)  was  calculated  initially  in 
order  to  obtain  sone  general  results  abar..  the  influence  of  the  deviation 
angle  on  the  flow  field.  For  this  case,  the  rotcr  operated  within  ar 
axisymmetric  boundary  layer  created  by  an  upstream  surface.  The  steps  ir. 
the  calculation  procedure  are  listed  in  Figure  50. 

The  measured  velocity  profile  in  the  rotor  plane  without  the  rotor 
was  extended  upstream  by  the  streamline  curvature  method  (Step  2)  in 
order  to  obtain  a velocity  profile  far  upstream  of  the  rotor.  This 
upstream  profile  is  needed  as  a boundary  condition  of  the  flow  field.  A 
comparison  between  the  measured  profile  in  the  rotor  plane  and  a profile 
obtained  by  projecting  the  upstream  profile  to  the  rotor  plane  using  the 
streamline  curvature  method  is  shown  in  Figure  56. 

The  theoretical  design  circulation  distribution  for  the  rotor  is 
shown  in  Figure  57.  Most  of  the  loading  occurs  near  the  mid-radius  where 
the  camber  angles  are  the  largest.  The  stagger  angles  are  high  as  a 
result  of  the  design  flow  coef f icier.t  being  less  than  one.  The  blade 
sections  are  thin  except  near  the  inner  wall.  The  blade  outlet  angles 
(S^)  were  determined  from  the  blade  geometry  by  using 

6 = >•  - X , (66) 

2 2 

where  \ is  the  stagger  angle  and  is  the  outlet  camber  angle.  A 
description  of  the  various  geometric  angles  is  shown  in  Figure  5b. 
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Figure  56.  Calculated  Meridional  Velocity  Profile 
for  Basic  Flow  No.  11 
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The  first  estimate  cf  the  flow  outlet  angles  (Step  3)  was  calculated 

using  the  blade  outlet  angles  (8  ) and  a deviation  angle  based  on 

2 

Howell's  correlation  given  by  Equation  (27).  The  flow  outlet  angles  are 
given  in  Table  5 for  Basic  Flow  No.  1. 

Using  the  estimated  flow  outlet  angles  and  the  design  flow  < inef- 
ficient, the  flow  field  was  calculated  using  the  streamline  curvature 
equations  (Step  4).  There  were  forty-eight  streamlines  calculated  in  the 
flow  field  of  which  twenty-eight  were  through  the  rotor  plane,  the  first 
streamline  being  at  the  inner  wall  and  the  twenty-eighth  streamline  going 
through  the  rotor  tip.  The  streamlines  through  the  rotor  were  spaced 
more  closely  near  the  inner  wall  because  secondary  flows  are  most  impor- 
tant near  the  inner  wall.  Some  of  the  results  at  several  radii  are  given 
in  Table  5.  The  streamline  curvature  equations  are  inviscid  so  that 
there  is  a finite  positive  velocity  at  the  inner  wall.  More  details  ol 
the  flow  field  for  all  the  streamlines  are  given  in  Reference  _/  55  /. 

The  velocity  components  given  .in  Table  5 are  for  radii  used  in  the 
flow  outlet  angle  calculations.  A comparison  shows  a large  acceleration 
of  the  axial  velocity  through  the  rotor  with  the  exception  of  the  rotor 
tip  where  the  flow  is  most  difficult  to  define. 

A second  estimate  of  the  flow  outlet  angles  (Step  5)  for  this  rotor 
was  made  using  the  flow  field  calculated  in  Step  4.  The  inlet  flow 
angles  (8  ) and  the  axial  velocity  ratios  (AVR)  were  computed  so  that  a 
method  developed  by  Lieblein  / 40_/  was  applied  to  estimate  the  deviation 

angle  for  any  blade  incidence  angle.  The  effects  of  axial  acceleration 

★ 

(A51),  blade,  thickness  (A<5  ),  and  camber  (<5  ) given  by  Equations  (29), 
(30),  and  (31),  respectively,  were  calculated  for  the  primary  blade  radii. 


a 
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TABLE  5 

OUTLET  ANGLES  USING  HOWELL'S  CORRELATION 
(STEP  4)  FOR  BASIC  FLOW  NO.  1 
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0.26 

0.401 

65.8 

0.489 

0.260 

52.3 

1.217 

0.30 

0.509 

63.7 

0.569 

0.254 

53.8 

1.117 

0.40 

0.631 

65.3 

0.719 

0.249 

57.4 

1.139 

0.50 

0.706 

67.7 

0.828 

0.229 

60.9 

1.172 

0.60 

0.764 

69.7 

0.861 

0.175 

65.5 

1.127 

0.70 

0.812 

71.4 

0.886 

0.129 

68.8 

1.091 

0.80 

0.841 

73.0 

0.895 

0.087 

71.4 

1.064 

0.90 

0.864 

74.4 

0.916 

0.055 

73.2 

1.061 

0.95 

0.872 

75.1 

0.854 

0.021 

75.2 

0.979 
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are  given  in  Table  6 for  the  primary  rotor  radii.  It  is  important  to 
note  that  the  major  differences  in  the  second  calculation  of  the  flow 
field  as  compared  to  the  first  calculation  with  Howell's  correlation  arc 
due  to  the  effects  of  acceleration  near  the  mid-radius  and  the  effects  of 
blade  thickness  near  the  inner  wall.  The  tangential  velocity  near  the 
inner  wall  was  reduced  primarily  by  the  effect  of  thickness  and  the  axial 
velocity  near  the  mid-radius  was  increased  by  including  the  AVK  effect. 

An  estimate  of  the  rotor  streamlines  was  calculated  from  the  flow 

i 

analysis  using  the  second  estimate  of  flow  outlet  angles.  The  effect  of 
secondary  flows  can  now  be  calculated  along  the  twenty-eight  rotor 

streamlines.  Equations  (37)  and  (41)  were  used  to  calculate  « , and  , 

2 *2 

respectively . 

[ 

These  secondary  flow  equations  give  the  change  in  vorticity  due  to 

the  presence  of  the  blades.  The  difference  in  the  streamwise  vorticity 

between  the  rotor  exit  (to  ,)  and  the  inlet  (w  ,)  is  the  passage  secondary 

2 1 

vorticity.  The  effect  of  this  additional  vorticity  was  found  by  initially 
solving  for  a secondary  stream  function  (Equation  55))  in  the  rotor  exit 
plane.  The  solution  of  Equation  (55)  was  then  used  in  Equation  (64)  to 
obtain  a deviation  angle  due  to  the  passage  secondary  vorticity  ( ’ ' ). 
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TABLE  6 

OUTLET  ANGLES  USING  NASA'S  CORRELATION 


(STEP  5) 

FOR 

BASIC  FLOW 

NO.  1 

e/rr 

B 

2 

6 

o 

* 

A6 

A6 ' 

* 

B 

2 

0.26 

49.0 

2.6 

10.1 

4.5 

57.2 

0.30 

49.8 

3.5 

6.9 

2.2 

58.0 

0.40 

52.7 

4 . 6 

4.2 

2.5 

59.0 

0.50 

55.9 

5.4 

2.5 

2.9 

60.9 

0.60 

59.6 

5.7 

1.7 

2.0 

65.0 

0.70 

63.5 

5.1 

1.3 

1.4 

68.5 

0.80 

67.2 

4.4 

0.7 

0.9 

71.4 

0.90 

70.2 

3.2 

0.6 

0.3 

73.2 

0.95 

73.3 

1.5 

0.3 

0.0 

75.1 

Inlet 


Outlet 
Velocity  Profile 


8/Rr 

V /V 

X «° 

6 

l 

V /V 
x 00 

V /V 

9 “> 

* 

6 

2 

AVR 

0.26 

0.370 

67.5 

0.432 

0.215 

57.5 

1.167 

0.30 

0.483 

64.9 

0.516 

0.209 

57.9 

1.069 

0.40 

0.618 

65.8 

0.692 

0.244 

58.9 

1.120 

0.50 

0.704 

67.8 

0.831 

0.229 

60.9 

1.181 

0.60 

0.768 

69.6 

0.877 

0.185 

64.9 

1.148 

0.70 

0.818 

71.3 

0.905 

0.139 

68.3 

1.107 

0.80 

0.845 

72.9 

0.899 

0.090 

71.3 

1.063 

0.90 

0.867 

74.4 

0.921 

0.058 

73.1 

1.063 

0.95 

0.874 

75.0 

0.860 

0.024 

75.1 

0.984 
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The  effect  of  the  passage  secondary  vorticitv  can  be  seen  in  Table 
7 which  lists  the  final  flow  outlet  angles  used  in  the  calculation  of 
the  final  flow  field  for  Basic  Flow  No.  1.  .Secondary  flows  are  only 
significant  near  the  inner  wall  where  the  inlet  velocity  gradients  are 
the  largest.  The  result  is  an  overturning  of  the  flow  field  near  the 
inner  wall  where  the  inlet  velocity  gradients  are  the  largest.  The 
result  is  an  overturning  of  the  flow  field  near  the  inner  wall  giving  a 
larger  tangential  velocity  (V  ) than  expected.  Outside  this  small  region 
near  the  inner  wall,  the  secondary  passage  vorticitv  does  not  have  a 
significant  effect  on  the  deviation  angle.  The  calculations  gave  a devia- 
tion angle  change  of  less  than  0.2°  in  this  outer  region. 

A third  calculation  of  the  flow  field  (Step  9)  was  done  by  including 
the  secondary  flow  effects  with  the  previous  NASA  deviation  angles. 

Because  each  deviation  angle  was  calculated  for  a slightly  different  ilcv 
field,  a check  was  made  on  the  effects  of  axial  acceleration,  camber,  and 
blade  thickness  on  the  deviation  angles  using  the  flow  field  calculated 
in  Step  9.  The  results  showed  changes  less  than  0.2°  in  the  final  devia- 
tion angle  (6).  Therefore,  the  flow  field  calculated  in  Step  9 was  the 
final  flow  field.  The  final  deviation  angles  for  Basic  Flow  No.  1 are 
given  in  Table  7. 

A correlation  can  now  be  made  between  the  calculated  flow  field 

(Step  9)  and  the  measured  rotor  inlet  and  outlet  profiles  which  were 

discussed  in  Section  3.3.  A comparison  between  the  calculated  and 

measured  inlet  velocity  profiles  to  the  rotor  is  shown  in  Figure  59.  In 

this  figure,  the  measured  meridional  velocity  profile  was  corrected  to 

the  axial  velocity  profile  using  the  streamline  angles  which  were  less 
o 

than  8 . This  figure  shows  quite  good  correlation  except  near  the 


1 
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inner  wall  where  the  inviscid  theory  over  predicts  the  axial  velocity. 

A comparison  between  the  calculated  and  measured  rotor  outlet  profile  is 

shown  in  Figure  60.  The  agreement  in  the  axial  velocity  profile  is  quite 

! 

good  including  the  effect  of  the  rotor <tip  shown  by  the  dip  in  the  axial 
velocity.  The  tangential  velocity  profiles  are  also  in  good  agreement. 

The  differences  near  the  rotor  tip  are  due  to  nonaxisymmetric  effects. 

Near  the  inner  wall,  the  calculated  profile  over  predicts  the  measured 
profile.  This  difference  can  be  either  attributed  to  the  inviscid  second- 
ary flow  theory  or  to  the  measurements  being  made  too  close  to  the  rotor 
exit  plane.  This  second  point  is  most  important.  The  secondary  flow 
theory  assumes  that  the  induced  velocity  field  occurs  in  the  rotor  exit 
plane  when,  in  fact,  the  secondary  vorticity  filaments  extend  downstream. 

5.1.2  Other  Basic  Flows.  The  calculation  procedure  shown  in  Figure 
50  for  determining  the  flow  field  for  a rotor  was  applied  to  other  basic 
flow  configurations.  Most  of  the  Basic  Flow  Nos.  given  in  Table  1 were 
calculated.  However,  the  results  for  Basic  Flow  Nos.  1,  2,  3,  and  4 are 
presented  because  complete  experimental  velocity  data  were  only  available 
for  these  flows. 

The  velocity  profiles  in  the  plane  of  the  rotor  without  the  rotor 
used  in  the  initial  calculation  of  each  basic  flow  were  obtained  from 
experimental  data  presented  in  Section  3.3.1.  All  of  the  basic  flows  had 
an  axisymmetric  velocity  profile  with  the  exception  of  Basic  Flow  No.  4. 

For  this  case  which  had  upstream  struts,  a circumferentially  averaged 
profile  was  used. 

* 

The  flow  outlet  angles  (B  ) were  determined  for  each  basic  flow  bv 

2 

die  method  previously  outlined.  The  final  angles  are  given  in  Table  7 
for  Basic  Flow  No.  1,  in  Table  8 for  Basic  Flow  No.  2,  in  Table  9 for 
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Figure  60.  Rotor  Outlet  Velocity  Comparison  for 
Basic  Flow  No.  1 
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TABLE  8 

FINAL  OUTLET  ANGLES  FOR  BASIC  FLOW  NO.  2 


R/Ra 

B 

O 

6 

0 

* 

A6 

A6 ' 

A 6 

s 

* 

8 

2 

0.26 

49.0 

2.6 

9.6 

3.4 

-5.9 

51.9 

0.30 

49.8 

3.4 

7.0 

1.1 

-3.2 

55.9 

0.40 

52.7 

4.7 

4.0 

2.2 

-1.1 

58.1 

0.50 

55.9 

5.4 

2.3 

2.9 

— 

60.7 

0.60 

59.6 

5.8 

1.7 

2.3 

— 

64.8 

0.70 

63.5 

5.2 

1.2 

1.7 

— 

68.2 

0.80 

67.2 

4.5 

0.8 

1.4 

— 

71.1 

0.90 

70.2 

3.2 

0.6 

1.1 

— 

72.9 

0.95 

73.3 

1.9 

0.3 

0.0 

— 

75.5 

Inlet 

Velocity  Profile 


r/rr 

V /V 
x * 

6 

1 

0.26 

0.465 

64.9 

0.30 

0.569 

63.6 

0.40 

0.693 

65.6 

0.50 

0.779 

67.8 

0.60 

0.843 

69.8 

0.70 

0.891 

71.6 

0.80 

0.912 

73.4 

0.90 

0.914 

75.1 

0.95 

0.913 

75.9 

Outlet 


Velocity  Profile 


V /V 

X °° 

Vv 

0 00 

* 

8 

2 

0.543 

0.301 

51.8 

0.594 

0.269 

55.9 

0.774 

0.280 

58.2 

0.917 

0.279 

60.7 

0.966 

0.231 

64.9 

0.998 

0.184 

68.2 

1.000 

0.132 

71.1 

1.034 

0.095 

72.8 

0.927 

0.043 

75.5 
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Basic  Flow  No.  3,  and  in  Table  10  for  Basic  Flow  No.  4.  A comparison 
between  basic  flows  shows  large  differences  in  outlet  angles  near  the 
inner  wall.  Additional  streamline  data  are  given  in  Reference  J_  55  /. 

A comparison  between  the  calculated  and  experimental  rotor  inlet  and 
outlet  profiles  is  given  in  Figures  59  and  60  for  Basic  Flow  No.  1,  in 
Figures  61  and  62  for  Basic  Flow  No.  2,  in  Figures  63  and  64  for  Basic 
Flow  No.  3,  and  in  Figures  65  and  66  for  Basic  Flow  No.  4.  In  general, 
the  calculations  are  in  good  agreement  with  the  experimental  results. 

The  calculated  axial  velocity  near  the  inner  W3ll  is  greater  for  the 
inlet  profiles  and  for  most  outlet  profiles.  The  tangential  velocity 
profile  is  over  predicted  for  all  basic  flows  near  the  inner  wall.  This 
is  due  to  the  secondary  vorticitv  model  which  assumes  that  the  influence 
of  the  vorticity  on  the  velocity  field  occurs  in  the  trailing  edge  plane 
of  the  rotor.  Comparisons  between  basic  flows  show  that  the  tangential 
velocity  increased  significantly  by  operating  at  a low’  flow  coefficient 
and  by  adding  upstream  struts. 

The  torque  coefficient  is  a measure  of  the  magnitude  of  the  tangen- 
tial velocity.  It  is  defined  as 


V 


Ri/RR 


'r^xf  “'t' 

OO  orj  R R 


(68) 


This  coefficient  was  calculated  from  the  flow  field  for  each  basic  flow 
and  was  normalized  by  the  coefficient  obtained  from  a reference  which  was 
Basic  Flow  No.  1. 


These  torque  results  obtained  from  the  calculated  flow  field  are 
compared  to  experimental  data  as  shown  in  Figure  67.  The  changes  in 
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FINAL  RESULTS  FOR  BASIC  FLOW  NO.  4 


k/Br 

3 

2 

5 

0 

* 

A6 

> r 1 
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A6 

s 

* 

2 

0.26 

49.0 

2.5 

8.7 
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-5.0 
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0.30 

49.8 

3.4 
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0.1 

-3.2 

56.5 

0.40 
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4.6 
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2.9 

-1.3 

57.2 

0.50 

55.9 

5.4 

2.5 

2.9 

— 

60.9 

0.60 

59.6 
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2.3 

— 

64.7 

0.70 
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1.6 

— 

68.3 

O.SO 
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0.8 
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— 
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— 
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0.95 
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— 
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x 00 

£ 

i 
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* 
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0.26 
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0.30 

0.542 

62.3 
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56.5 

0.40 

0.608 

66.1 

0.709 

0.270 

57.3 

0.50 

0.704 

67.8 

0.1319 

0.251 

60.9 

0.60 

0.759 

69.8 

0.373 

0.209 

64.8 

0.70 

0.802 

71.6 

0.896 

0.162 

68.3 

0.80 

0.827 

73.3 

0.383 

0.111 

71.5 

0.90 

0.834 

74.9 

0.952 

0.091 

72. 4 

0.95 
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0.835 
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Figure  63.  Rotor  Inlet  Velocity  Comparison  for 
Basic  Flow  No.  3 


EXPERIMENTS 

DATA 


1.0  1.2 


VELOCITY  RATIO,  V /V 

X °o 


Figure  65.  Rotor  Inlet  Velocity  Comparison  for 
Basic  Flow  No.  4 
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torque  due  to  adding  struts  or  operating  off-design  are  predicted  very 
well . 

5.1.3  Comparisons  Between  Calculated  and  Measured  Flow  Fie 1 ds  at 
the  End  of  the  Rotor  Cap.  So  far,  only  comparisons  between  flow  field 
calculations  and  flow  measurements  near  the  rotor  plane  have  been  pre- 
sented. The  measured  tangential  velocity  profile  near  the  rotor  plane 
does  not  show  a well  defined  vortex  near  the  inner  wall.  However, 
measurements  of  the  tangential  velocity  profile  at  the  end  of  the  cap  do 
show  a vortex. 

Figure  68  shows  a comparison  at  the  end  of  the  rotor  cap  between  the 
final  flow  field  calculation  (Step  9)  and  experimental  data  for  Basic 
Flow  No.  1.  The  tangential  velocity  profiles  are  in  very  good  agreement 
outside  of  the  vortex  core:  however,  the  axial  velocity  profiles  do  not 
agree  near  the  vortex  core.  This  is  opposite  to  the  correlation  trends 
obtained  at  the  rotor  outlet  between  the  calculated  flow  field  and 
experimental  data. 

The  end  of  the  rotor  cap  is  located  several  chord  lengths  downstream 
of  the  rotor  exit  plane.  Therefore,  it  is  not  surprising  that  the  agree- 
ment between  the  calculated  flow  field  including  secondary  flows  and 
measured  tangential  velocity  profile  agreement  is  better  at  the  end  of 
the  rotor  cap.  The  tangential  velocity  profile  calculated  by  the  stream- 
line curvature  equations  using  a secondary  flow  deviation  angle  compares 
well  with  the  tangential  velocity  profile  data  outside  of  the  rotational 
region  of  the  vortex.  A poor  agreement  is  found  inside  the  vortex  core 
because  of  the  assumptions  used  in  the  development  of  the  streamline 
curvature  equation. 


Tn  addition  to  the  influence  of  the  passage  streamwise  secondary 
vorticity  on  the  tangential  velocity  profile,  the  secondary  voriicity 
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also  changes  the  axial  velocity  profile.  However,  the  net  mass  flow- 
change  to  the  axial  velocity  profile  due  to  secondary  vorticity  must  be 
zero . 

The  calculation  of  the  change  in  axial  velocity  occurring  downstream 
of  a rotor  plane  due  to  the  secondary  vorticity  is  a difficult  problem. 
Even  assuming  that  downstream  of  the  blade  row  the  vorticity  become 
distributed  circumferentially,  few  good  correlations  with  axial  velocity 
exists  and,  in  most  cases,  axial  velocity  changes  are  ignored. 

Hawthorne  _/  42 _/  has  shown  that  there  are  three  components  of  axial 
vorticity  downstream  of  the  blade  exit  plane.  One  is  the  passage  sec- 
ondary vorticity;  the  other  two  are  the  trailing  shed  vorticity  and  the 
trailing  filament  vorticity  which  lie  along  the  blade  wake.  The  trailing 
shed  vorticity  is  caused  by  a gradient  in  circulation  along  the  blade 
span  and  the  contributions  of  both  the  primary  and  secondary  parts  are 
included  in  the  final  flow  field  analysis  (Step  9).  The  secondary  part 
of  the  trailing  shed  vorticity  is  the  change  in  blade  turning  due  to  the 
passage  secondary  vorticity. 

Trailing  filament  vorticity  arises  from  the  cellular  motion  induced 
in  the  blade  passage  by  the  secondary  streamwise  vorticity.  Dixon  J_  50 _/ 
derives  an  expression  for  the  axial  component  of  the  combined  passage  and 
trailing  filament  vorticities.  The  results  are 

* N 

w = (<d  , - u>  ,)cosf<  - u : — , (f>9) 

x.  s'  s'  ? o irr 

A 2 1 


where  u is  the  radial  perturbation  velocity  due  to  the  secondary  passage 
o 

vorticity  evaluated  at  0 = 0,  2tt/N,  ...  . The  radial  velocity  is  given 
by  Equation  (49),  knowing  the  solution  of  the.  stream  function  given  by 
Equation  (55). 

The  axial  component  of  vorticity  can  be  integrated  over  the  radial 
direction  to  give  an  axial  velocity  variation.  The  constant  of  integra- 
tion is  obtained  from  continuity.  Again,  the  net  mass  flow  through  a 
downstream  plane  induced  by  the  axial  component  of  vorticity  must  be  zero. 

Calculation  of  <i»x  for  Basic  Flow  Ho.  1 using  Equation  (69)  shows  that 
the  axial  secondary  vorticity  is  approximately  zero  over  the  blade  span. 
Thus,  the  axial  component  of  passage  secondary  vorticity  and  the  axial 
component  of  trailing  filament  vorticity  are  approximately  equal  and 
opposite  at  each  streamline.  Therefore,  for  Basic  Flow  No.  1 where  the 
rotor  operates  within  the  inlet  velocity  gradient,  no  variation  in  axial 
velocity  profile  due  to  axial  vorticity  would  be  anticipated.  This  result 
was  also  shown  by  Smith  j_  56_/. 

However,  the  aforementioned  result  does  not  explain  the  large  dif- 
ferences between  measured  and  calculated  axial  velocity  near  the  vortex 
core.  The  measured  axial  velocity  profile  correlated  well  with  measure- 
ments made  near  the  rotor  plane.  One  possible  hypothesis  for  this 
discrepancy  is  that  the  secondary  flow  field  does  not  become  distributed 
downstream  of  the  rotor  plane.  Although  the  net  axial  secondary  vorticity 
remains  approximately  zero,  the  passage  vorticity  remains  concentrated  and 
dominates  the  flow  field  near  the  rotor  inner  wall. 

For  this  flow  field,  the  axial  velocity  variation  calculated  by 
Equations  (47)  and  (48)  used  in  the  solution  of  the  secondary  stream 
function  was  applied  to  the  calculated  profile  at  the  end  of  the  rotor 
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cap.  The  variation  calculated  for  each  of  the  twenty-eight  streamlines 
in  the  rotor  exit  plane  was  added  to  the  corresponding  streamline  in  the 
rotor  cap  plane.  The  results  are  shown  in  Figure  69  for  Basic  Flow  No.  ]. 
The  correlation  between  measured  and  calculated  axial  velocity  profiles  is 
very  good  using  only  the  passage  secondary  vorticity. 

Additional  comparisons  between  the  final  flow  field  calculations 
(Step  9)  and  experimental  data  obtained  at  the  rotor  cap  are  given  in 
Figures  70,  71,  and  72.  Comparisons  between  measured  and  calculated 
velocity  profiles  for  Basic  Flows  No.  2,  3,  and  4 are  shown  in  Figures  70, 
71,  and  72,  respectively.  In  all  of  these  figures,  an  axial  velocity 
variation  due  to  secondary  axial  vorticity  was  not  included. 

5 . 2 Calculation  of  Secondary  Vorticity  and  Vortex  Parameters 

5.2.1  Secondary  Vorticity  for  Basic  Flow  No. 1.  In  the  schematic 

for  the  calculation  of  the  final  flow  field  shown  in  Figure  50,  the  secon- 
dary flow  calculation  is  shown  in  Step  7.  This  calculation  not  only  adds 
the  blade  to  blade  flow  effects  due  to  secondary  vorticity  to  the  final 
flow  field  but  also  gives  the  vortex  parameters.  The  substeps  of  Step  7 
are  shown  in  Figure  73. 

For  Step  7B,  the  vorticity  was  numerically  calculated  along  a rela- 
tive streamline  (s')  using  Equation  (34)  for  the  relative  normal  component 
of  the  absolute  vorticity  (w^,)  and  Equation  (41)  for  the  relative  stream- 
wise  component  of  the  absolute  vorticity  (w  , ) . The  flow  parameters 
necessary  for  the  vorticity  calculation  were  defined  along  a relative 
streamline  from  the  flow  field  calculation  in  Step  6.  In  all,  the 
vorticity  for  twenty-eight  relative  streamli;  s was  calculated  for  the 

rotor.  The  initial  vorticity  components  (u>  ,,  w .)  were  calculated  from 

s n 

I 1 
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STEP  7 SECONDARY  FLOW  CALCULATION 


STEP  7A 


STEP  7B 


STEP  7C 


STEP  7D 


STEP  7E 


Figure  73. 


Outline  of  Secondary  Flow  Calculations 


the  inlet  velocity  gradient  to  the  rotor  as  shown  in  Figure  54.  The 
rotor  inlet  flow  was  axisymmetric  and  the  velocity  was  nearly  axial. 
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Therefore,  the  normal  component  of  absolute  vorticity  (w  ) is  the  gradient 

■ n 

of  the  meridional  velocity  profile.  The  normalized  normal  component  of 

vorticity  (w  ) is  shown  in  Figure  6 for  Basic  Flow  No.  1. 
n 

1 

As  shown  in  Equation  (34),  the  relative  component  of  the  normal 
vorticity  depends  on  the  variation  of  relative  velocity  (W)  along  a rela- 
tive streamline,  on  the  relative  component  of  the  absolute  inlet  normal 

vorticity  (w  ,),  and  on  the  flow  convergence-d i vergence  in  the  bi -normal 
1 

direction  (a,  ,).  The  results  for  Basic  Flow  No.  1.  are  given  in  Figure  74. 
b 

The  exit  normal  vorticitv  (w  ,)  is  not  very  different  from  the  inlet 

n 

2 

normal  vorticity  (w  ,). 

1 

The  normal  vorticity  calculations  were  used  in  Equation  (41)  for  the 
calculation  of  the  relative  streamwise  vorticity.  The  results  for  Basic 
Flow  No.  1 are  given  in  Figure  75.  The  exit  relative  streamwise  vorticity 
(u  ,)  is  very  different  from  the  inlet  relative  streamwise  vorticitv 

<v>- 

1 

The  terms  of  Equation  (41)  are  shown  in  Figure  76.  The  first  term  is 
simply  the  turning  of  a normal  component  of  vorticity  and  is  small  when 
compared  to  the  other  terms.  The  second  term  is  also  small  and  is  the 
curvature  induced  vorticity.  The  third  term  is  large  and  is  due  to 
rotation-induced  vorticity.  This  result  can  be  anticipated  because  the 
rotor  operates  at  a low  flow  coefficient  and  hence  is  large.  The  last 
term  is  simply  the  stretching  of  the  vorticity  by  the  flow  and  is  large 
because  of  the  high  stagger  angles  at  which  the  rotor  operates. 

The  secondary  passage  vorticity  is  the  difference  between  the  inlet 
and  outlet  vorticities  and  is  used  in  Step  7C  and  7D  for  the  calculation 


WITHOUT  UPSTREAM  STRUTS 
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DESIGN  FLOW  COEFFICIENT 
(BASIC  FLOW  NO.  1) 


Figure  76.  Components  of  the  Relative  Streamwise  Vortieity  for 
Basic  Flow  No.  1 
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of  the  deviation  angle  (A5g)  due  to  secondary  flows.  The  vorticity 

results  are  shown  in  Figure  77  for  Basic  Flow  No.  1.  The  normal  secondary 

vorticitv  (u  ,-w  ,)  is  very  small  and  is  neglected.  This  result  has  been 
n n 
2 1 

shown  by  many  secondary  flow  investigations.  The  streamwise  secondary 
vorticity  is  large  and  positive.  This  positive  vorticity  causes  an 
overturning  of  the  flow  near  the  inner  wall  for  this  rotor.  The  vorticity 
is  concentrated  near  the  inner  wall  which  correlates  with  the  large  devia- 
tion angles  calculated  near  the  inner  wall  and  the  very  small  deviation 
angles  calculated  over  the  remainder  of  the  blade  span. 

The  most  interesting  result  from  the  calculation  of  the  streamwise 

vorticity  is  the  shape  of  the  exit  vorticity  curve  (uj  ,)  shown  in  Figure 

2 

75.  There  is  a large  amount  of  negative  vorticity  near  the  inner  wall 
and  very  little  vorticity  in  the  mid-radius  region  of  the  rotor. 

5.2.2  Calculation  of  Vortex  Parameters.  The  passage  streamwise 
vorticity  (Step  7B)  was  calculated  from  Basic  Flow  Nos.  1,  2,  3,  and  4. 

The  shapes  of  the  curves  (u  ,)  were  all  similar  except  that  the  region  of 

2 

negative  vorticity  near  the  inner  wall  was  different  in  size  and 
magnitude . 

This  vorticity  near  the  inner  wall  can  coalesce  in  the  blade  passage 

to  form  a vortex-type  flow  as  shown  by  Reference  ]_  5 7 / . In  this  case, 

the  vorticity  downstream  of  the  rotor  will  be  nonaxisymmetric . However, 
the  vorticity  can  also  remain  distributed  in  the  passage  to  form  down- 
stream a vortex  with  axisymmetric  vorticity.  In  either  case,  there  will 

be  a characteristic  length  (L  ) and  circulation  (f)  associated  with  the 

c 

passage  vorticity  which  controls  the  minimum  pressure  of  the  resulting 
vortex  system.  The  vorticity  will  still  induce  a tangential  velocity 
distribution  similar  to  that  measured  cn  the  rotor  cap.  However,  large 
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fluctuations  in  the  tangential  velocity  would  be  expected  if  the  vorticlty 
remains  nonaxisymmetr ic . 

As  can  be  seen  in  Figure  75,  the  streamwise  vorticitv  (ui  ,)  near  the 

2 

inner  wall  has  a characteristic  length  (L  ) associated  with  it.  This 

c 

length  is  defined  as  the  difference  between  the  radius  where  the  vortlcity 
changes  sign  and  the  inner  wall  radius. 

A measure  of  the  circulation  associated  with  this  vorticity  was  found 
by  integrating  the  vorticity  (Step  7E)  from  the  inner  wall  to  the  radius 
where  the  vorticity  changes  sign.  The  equation  is 


where  rcosB  is  normal  to  the  vorticity.  The  calculated  circulation  is 

2 

per  blade  passage. 

The  vortex  parameters  (Step  7E)  for  the  four  basic  flow  configura- 
tions are  given  in  Table  11.  Also,  the  nondimensional  ratio  (F/L  V )) 

C 00 


which  is  a measure  of  the  minimum  pressure  coefficient  of  a Rankine  vortex 


is  given.  In  addition,  the  planar  momentum  thickness  of  the  mean  velocity 
profile  entering  the  rotor  is  given  as  an  indication  of  the  velocity 
deficit . 


5. 3 Correlation  of  Vortex  Parameters  with  Cavitation  Data 

Becuase  of  the  uncertainty  of  the  final  vortex  structure  at  the  end 

of  the  rotor  cap,  an  absolute  calculation  of  C in  this  region  would  be 

P • 
rmin 

very  difficult.  However,  the  vorticity  calculated  in  the  blade  passage 
eventually  organizes  the  low  momentum  fluid  near  the  inner  wall  into  the 


final  vortex  structure.  In  addition,  the  vortex  exists  near  the  inner 


wall  where  the  pressure  in  this  region  is  not  only  controlled  by  the 
local  vortex  but  also  by  the  rotor  velocity  field  outside  of  the  vortex. 

The  vorticity  does,  however,  dominate  the  flow  field  near  the  inner 
wall  when  a rotor  is  designed  to  be  unloaded  in  this  region.  In  this 
case,  changes  in  this  vorticity  will  directly  influence  the  resulting 
minimum  pressure  of  the  vortex.  Therefore,  Equation  (10)  can  be  written 


ff  r ) 

pmin . j L V j 


A = P-  c ^ A 

C ; r _ i 

Pmin  j I 

B l.  l.LcV«>-l 


where  T and  are  the  characteristic  vortex  parameters  of  the  passage 
streamwise  vorticity.  The  letters  A and  B refer  to  flow  states.  Also, 
it  is  assumed  that 


R “L 
c c 


In  order  to  relate  the  C of  the  vortex  to  the  cavitation  number, 

pmin 

several  criteria  must  apply.  First  of  all,  the  cavitation  data  must  be 
free  of  gas  effects  as  discussed  in  Section  3.4.3.  This  is  to  insure 


a = -C 
d p . 

rmin 


Secondly,  the  minimum  pressure  in  the  final  vortex  structure  must  be 
controlled  by  the  streamwise  vorticity  associated  with  the  vortex  and  not 
by  the  pressure  drop  in  the  radial  plane  created  by  the  primary  flow 
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outside  of  the  vortex.  Also,  variations  in  inlet  velocity  gradient  must 
be  small  enough  not  to  significantly  affect  the  mean  rotor  loading.  All 
of  these  criteria  seem  very  restrictive;  however,  performance  of  most 
rotors  operating  within  a velocity  gradient  meets  these  requirements. 

A correlation  between  calculated  minimum  pressure  coefficients  and 

cavitation  data  can  be  made  using  Equation  (71)  given  a reference  point 

and  knowing  the  effect  of  Reynolds  number  on  the  passage  vorticity.  A 

reference  point  for  Basic  Flow  No.  1 of  o = 2.8  at  a velocity  of  15  ft/sec 

was  chosen.  The  influence  of  Reynolds  number  on  the  passage  vorticity 

(w  ,)  had  to  be  solved, 
s 

2 

As  noted  in  Table  1,  Basic  Flow  Nos.  .14,  15,  16,  and  17  are  special 
theoretical  cases  which  were  used  in  the  calculation  of  Reynolds  number 

effects  on  the  passage  vorticity  (co  , ) . For  Basic  Flow  No.  14,  the 

2 

calculation  procedure  given  in  Figure  50  was  followed.  A measured 
velocity  profile  in  the  rotor  plane  without  the  rotor  obtained  at  a 
Reynolds  number  of  3 x 10^  was  used  in  Step  2 of  the  calculation.  For 
Basic  Flow  No.  15,  the  calculation  procedure  was  again  followed,  but  using 
a measured  velocity  profile  in  the  rotor  plane  without  the  rotor  obtained 
at  a Reynolds  number  of  8.8  x 105  in  Step  2.  The  two  velocity  profiles 
are  shown  in  Figure  78. 

One  point  should  be  made  of  these  additional  calculations.  The 
normal  component  of  absolute  vorticity  (u^)  used  in  the  vorticity  calcula- 
tions for  Basic  Flow  Nos.  14  and  15  was  obtained  by  taking  the  derivative 
of  the  theoretical  inlet  velocity  profile  from  Step  6.  As  can  be  noted  in 
Figure  59,  the  measured  inlet  velocity  profile  does  not  always  agree  with 
the  theoretical  calculations.  Therefore,  Basic  Flow  No.  1 is  different 
from  Basic  Flow  No.  14  because  the  normal  vorticity  for  Basic  Flow  No.  1 


L 
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Figure  78.  Axial  Velocity  Profiles  for  Different 
Reynolds  Numbers 
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was  obtained  from  a measured  velocity  profile  in  front  of  the  rotor. 

Both  basic  flows  used  the  same  measured  velocity  profile  in  the  rotor 
plane  without  the  rotor  for  Step  2.  The  resulting  vortex  parameters  of 
r and  L were  sLightly  different  for  Basic  Flow  No.  1 as  compared  to 
Basic  Flow  No.  14. 

Basic  Flow  Nos.  16  and  17  were  calculated  in  order  to  show  that  the 

vorticity  calculation  based  on  a theoretical  velocity  profile  will  not 

influence  the  ratios  of  vortex  parameters  obtained  from  Equation  (71). 

Basic  Flow  No.  16  is  similar  to  Basic  Flow  No.  2.  Using  the  vortex 

parameters  in  Equation  (71),  the  C r..tio  of  Basic  Flow  No.  2 to  Basic 

n . 

' min 

Flow  No.  1 was  found  to  be  within  2%  of  the  C ratio  of  Basic  Flow  No. 

p . 

rmin 

16  to  Basic  Flow  No.  14.  In  all  cases,  the  vortex  parameters  of  F and 

were  smaller  using  the  theoretical  normal  vorticity  (w  ) as  compared  to 

1 

using  normal  vorticity  (m  ) obtained  from  experimental  data. 

1 

Now  using  Equation  (71)  with  a reference  point  for  Basic  Flow  No.  1, 

comparisons  between  cavitation  data  and  a calculated  C were  made. 

P • 
min 

Some  of  the  correlations  are  shown  in  Figures  79,  80,  and  81. 

Referring  to  Figure  79,  a correlation  of  Basic  Flow  No.  1,  Basic  Flow 
No.  4,  and  Reynolds  number  with  cavitation  data  is  given.  For  this  cor- 
relation, the  ratio  of 


min  j. 


#15 


= 1.30 


; rain. 


#14 


(73) 


was  used  to  calculate  a point  of  a = 3.64  for  a velocity  of  45  ft/sec.  A 
line  was  drawn  between  these  two  points  for  the  correlation  of  data 


obtained  for  Basic  Flow  No.  1.  The  ratio  of 


WITHOUT  SCREEN 
DESIGN  FLOW  COEFFICIENT 

o WITHOUT  UPSTREAM  STRUTS  (BASIC  FLOW  NO.  1) 
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Figure  79.  Correlation  with  Cavitation  Data  for 
Basic  Flow  Nos.  1 and  A 
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Equation  (75)  gives  the  effect  of  adding  upstream  screen.  Again,  the 
correlation  between  cavitation  data  and  the  calculated  vortex  parameters 
is  quite  good. 

A final  correlation  between  calculated  data  and  cavitation  data  is 
shovm  in  Figure  81.  This  figure  is  a composite  of  the  four  basic  flow 
configurations.  The  calculated  effect  of  struts,  a screen,  and  off-desigr. 
flow  coefficient  was  combined  to  give  the  calculated  curves.  Good  agree- 
ment was  found  except  for  the  effect  of  off-design  flow  coefficient. 

As  can  be  seen  in  Figures  42,  45,  and  46,  the  cavitation  data  shows 
a rapid  change  in  cavitation  number  with  velocity  for  off-design  flow 


Figure  80.  Correlation  with  Cavitation  Data  for 
Basic  Flow  Nos.  1,  3,  4,  and  5 
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coefficient.  Secondary  flow  theory  alone  cannot  account  for  this  trend. 
Therefore,  these  trends  must  include  some  additional  primary  flow  effects. 

A comparison  can  be  now  made  between  the  vortex  parameters  (T,  R ) 
obtained  from  flow  measurements  at  the  end  of  the  rotor  cap  and  the  vortex 
parameters  calculated  from  the  passage  exit  vorticity.  As  shown  in  Figure 
30,  solid  curves  can  be  drawn  through  the  tangential  velocity  distribution 
(V  /V^)  outside  of  the  viscous  region  and  inside  the  viscous  region.  The 
intercept  of  these  two  solid  curves  is  the  approximate  radius  of  the 
vortex  cone.  The  circulation  associated  with  this  core  region  is 

T = ttR  V , (76) 

c c 

where  is  the  tangential  velocity  at  the  intercept  of  the  two  curves. 

The  resulting  vortex  parameters  (T,  R^)  are  given  in  Table  12. 

Also  shown  in  Table  12  are  two  sets  of  vortex  parameters  (f,  L ) 
calculated  from  the  exit  passage  vorticity.  The  first  set  given  in 

Column  2 is  the  parameters  obtained  from  the  u , calculations  which  were 

2 

used  in  the  correlation  of  cavitation  data.  The  second  set  given  in 
Column  3 was  obtained  by  integrating  the  axial  passage  vorticity  given  as 

* * 

e)  = uj  , cos  8 + w . sin  3 (77) 

x S 2 n 2 

2 2 

over  the  characteristic  length  (L  ) defined  from  the  u , calculations. 

c s 

O 

A comparison  shows  a good  correlation  between  the  parameters  estima- 
ted from  the  data  and  the  parameters  calculated  in  Column  3.  In  par- 
ticular, the  axial  circulation  in  the  measured  core  region  appears  to  be 

the  circulation  obtained  within  the  characteristic  length  (L  ) defined  at 

c 


r 
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the  rotor  exit  plane.  In  addition,  the  core  of  the  vortex  measured  at 
the  end  of  the  cap  is  approximately  one-half  of  the  characterist ic  length 


6.  SUMMARY  AND  CONCLUSIONS 


The  critical  cavitation  number  of  a vortex  formed  along  the  inner 
wall  of  a rotor  operating  within  a boundary  layer  has  been  correlated 
with  the  streamwise  vorticity  calculated  in  the  blade  passage.  A method 
has  been  outlined  which  can  be  used  to  calculate  the  passage  streamwise 
vorticity  for  a lightly  loaded  rotor  having  a low  hub-to-tip  ratio.  The 
method  involves  a solution  to  the  secondary  flow  equations  along  a mean 
streamline  in  the  rotor  passage.  The  streamlines  were  estimated  using 
streamline  curvature  method  to  solve  the  equations  of  motion  in  the 
direct  turbomachinery  problem. 

An  outline  of  the  method  used  to  calculate  the  flow  field  including 
the  secondary  flow  calculations  is  given  in  Figure  50.  The  flow  field 
calculation  (Step  6)  using  the  streamline  curvature  method  gives  the  data 
necessary  for  the  calculation  of  the  secondary  flows.  In  particular, 
these  data  include  the  inlet  velocity  profile,  the  outlet  velocity 
profile,  and  the  streamline  positions.  The  normal  component  of  absolute 
vorticity  (uO  is  calculated  from  the  inlet  velocity  profile  using 
Equation  (16).  These  data  are  in  an  absolute  coordinate  system. 

The  secondary  flow  calculations  are  outlined  in  more  detail  in 
Figure  73.  A summary  of  the  equations  used  in  the  secondary  flow  calcu- 
lations is  as  follows: 


(1)  The  initial  flow  field  data  are  changed  into  a rotating 
coordinate  system.  Now,  the  relative  velocity  (W)  and 
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the  streamline  angle  relative  to  the  axial  direction  (6) 
are  known  at  any  position  along  the  relative  streamline. 

(2)  As  shown  in  Figure  54,  the  initial  relative  components 

of  absolute  vorticity  (u>  ,,  u ,)  are  calculated.  In 

s n 
1 1 

addition,  the  relative  components  of  the  absolute  rota- 
tion (S2  ,,  H ,,  12.  , ) are  calculated  along  the  relative 
s n b 

streamlines,  the  radius  of  curvature  of  the  relative 

■k 

streamline  (R  ) is  calculated  along  the  relative  stream- 
lines using  Equation  (43),  and  the  relative  bi-normal 
component  of  absolute  vorticity  (oj^,)  is  calculated  along 
a relative  streamline  using  Equation  (44) . 

(3)  The  relative  normal  component  of  absolute  vorticity 
(u^,)  is  calculated  along  each  relative  streamline  using 
Equation  (37). 

(4)  The  relative  streamwise  component  of  absolute  vorticity 
(ug))  is  calculated  along  each  relative  streamline  using 
Equation  (41).  The  secondary  vorticity  at  the  rotor 

exit  plane  for  each  streamline  Is  u ,-w  ,. 

2 1 

(5)  The  boundary  conditions  for  the  secondary  stream  function 
defined  by  Equation  (55)  are  calculated.  First,  the 

, secondary  velocity  perturbation  in  the  axial  direction 

(w)  is  calculated  using  Equations  (47)  and  (48). 

Secondly,  the  right-hand  side  of  Equation  (55)  is 


calculated . 
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(6)  A series  solution  for  the  secondary  stream  function  is 
obtained  for  Equation  (59)  after  the  calculation  of 
4^(z)  given  by  Equation  (60). 

(7)  The  secondary  velocity  perturbation  in  the  tangent i .1 
direction  (v)  is  calculated  using  Equation  (49)  and 
Equation  (63)  . 

(8)  The  deviation  angle  due  to  the  secondary  vorticity 

(AS  ) is  calculated  using  Equation  (65) . 
s 

The  results  of  the  secondary  flow  calculations  showed  that  the  exit 

streamwise  component  of  vorticity  (to  ,)  for  the  rotor  depended  upon 

2 

rotation  induced  secondary  vorticity.  In  addition,  the  direction  of  the 

passage  secondary  vorticity  (w  ,-w  ,)  was  opposite  to  that  of  the  exit 

2 1 

streamwise  vorticity  (w  _ , ) . The  passage  secondary  vorticity  for  the  rotor 

2 

caused  an  overturning  of  the  flow  near  the  inner  wall. 

Extensive  experimental  data  were  obtained  not  only  to  verify  the 

calculation  procedure  but  also  to  help  in  the  analysis  of  the  flow  field. 

Some  conclusions  can  be  made  from  the  experimental  daca  about  the  flow 

field  near  the  rotor  and  the  vortex  structure. 

For  the  flow  field,  these  are: 

(1)  As  the  inlet  velocity  profile  near  the  inner  wall  was 
increased,  the  cavitation  number  also  increased. 

(2)  The  momentum  thickness  of  the  inlet  velocity  profile  does 
not  directly  correlate  with  the  core  size  of  the  vortex. 

(3)  Changes  of  less  than  12%  in  inlet  velocity  profile  near 
the  inner  wall  caused  changes  as  much  as  61%  in  cavita- 


tion number  of  the  vortex. 
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(4)  Tangential  velocity  measurements  at  the  rotor  cap  show  a 
vortex  within  a constant  tangential  velocity  field  due 
to  the  rotor. 

The  conclusions  related  to  the  vortex  structure  are: 

(1)  The  maximum  tangential  velocity  measured  at  the  vertex 
core  radius  does  not  correlate  with  critical  cavitation 
data  of  the  vortex  assuming  chat  all  of  the  vor  deity  is 
concentrated  in  the  vortex  core. 

(2)  The  minimum  pressure  of  the  vortex  system  does  not  occur 
on  the  surface  of  the  rotor  cap  above  a critical  Reynolds 
number . 

(3)  Changes  of  inlet  velocity  profile  of  less  than  12%  cause 
the  vortex  core  radius  to  change  by  as  much  as  42%. 

(4)  The  vortex  core  size  was  found  to  be  inversely  propor- 
tional to  the  cavitation  number  of  the  vortex. 

Additional  conclusions  can  be  made  about  secondary  flow  generated 
vortex  cavitation  from  comparisons  of  flow  field  calculations  with  flow 
field  measurements.  These  are: 

(1)  The  vortex  parameters  obtained  from  the  calculated  exit 
streamwise  vorticity  profile  do  correlate  well  with 
critical  cavitation  data. 

(2)  Even  though  the  theory  was  inviscid,  the  correlations 
between  the  calculated  final  flow  field  and  measurements 
are  quite  good,  including  those  obtained  at  the  rotor 
cap . 

(3)  The  critical  cavitation  number  of  the  vortex  increased 
with  the  addition  of  upstream  struts,  with  increasing 
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Reynolds  number,  and  with  a 10%  reduction  in  coefficient 
below  the  design  value. 

(4)  The  critical  cavitation  number  of  the  vortex  decreased 
with  the  addition  of  a screen  and  with  a 10%  increase  in 
flow  coefficient  above  the  design  value. 

(5)  The  effects  of  Reynolds  number,  struts,  screens,  and 
off-design  flow  coefficients  can  be  combined  to  success- 
fully correlate  cavitation  data  for  any  basic  flow 
considered . 

(6)  The  secondary  flow  theory  with  its  many  assumptions  gives 
a good  prediction  of  the  tangential  velocity  profile  near 
the  inner  wall. 

No  single  calculation  or  single  piece  of  experimental  data  can 
explain  the  vortex  structure.  However,  collectively,  a model  of  the  vor- 
tex system  becomes  apparent  above  a critical  Reynolds  number.  The  fact 
that  the  calculated  passage  vorticity  and  its  characteristic  length  scale 
correlates  the  cavitation  data  and  that  the  measured  vortex  core  radius 
and  maximum  tangential  velocity  do  not  correlate  the  cavitation  data 
indicates  that  the  net  circulation  of  the  vortex  is  not  within  the  vortex 
core.  Thus,  the  complicated  flow  around  the  rotor  cap  appears  to  be  like 
that  of  several  small  vortices  inside  of  larger  single  vortex  where 
viscous  effects  do  not  dominate  the  flow  near  the  inner  wall  of  the  rotor. 
Each  one  of  the  smaller  vortices  is  formed  within  each  blade  passage. 

This  flow  picture  explains  the  vortex  shown  in  Figure  1. 

This  investigation  of  secondary  flow  generated  cavitation  is  only  a 


beginning.  More  experimental  data  are  needed  of  the  secondary  flow  field 
in  the  rotor  passage  to  establish  the  motion  predicted  by  the  theoretical 
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development.  In  addition,  the  vorticity  roll-up  phenomenon  must  be 
investigated  in  order  to  predict  the  true  vortex  core  size  before  an 

absolute  calculation  of  C can  be  made.  Future  plans  are  to  lncor- 

P j 
rmin 

porate  the  secondary  flow  theory  into  the  design  of  a rotor  in  order  to 
obtain  a vortex  of  specified  critical  cavitation  index. 
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APPENDIX  A:  ANALYSIS  OF  EXPERIMENTAL  ERROR 


In  the  experimental  part  of  this  vortex  study,  many  different 
instruments  were  used.  The  following  is  a summary  of  the  measurements 
that  were  made  and  the  instruments  that  were  used: 

48-Inch  Wind  Tunnel 

(L)  Rotor  inlet  velocity  profiles  were  calculated  from  pressure 
data  obtained  from  separate  total  and  static  pressure 
probes. 

(2)  Static  pressure  taps  on  the  surface  were  used  to  calculate 
the  pressure  distribution  on  the  stationary  surface  up- 
stream of  the  rotor  plane. 

48-Inch  Water  Tunnel 

(1)  Rotor  outlet  velocity  profiles  were  calculated  from  data 
obtained  with  a 5-hole  prism  probe. 

(2)  Static  pressure  taps  on  both  the  stationary  surface  and 
rotating  core  were  used  to  calculate  the  surface  pressure 
distribution. 

(3)  Desinent  cavitation  numbers  for  the  vortex  were  calculated 
from  pressures  obtained  from  a pitot-static  probe  located 
in  the  free-stream. 

(4)  Velocity  profiles  were  calculated  at  the  end  of  the  rotor 
cap  from  frequency  data  of  a LDA. 

(5)  Several  rotor  inlet  and  outlet  velocity  profiles  were 
calculated  from  LDA  data. 

The  determination  of  accuracy  is  difficult  when,  as  in  most  cases, 
only  one  instrument  was  used  with  a limited  number  of  samples.  However, 
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the  rotor  outlet  tangential  velocity  profiles  for  several  cases  were 
measured  with  a 5-hole  prism  probe  and  the  LDA.  These  two  instruments 
were  used  during  different  water  tunnel  test  programs.  The  results  show 
an  experimental  accuracy  of  3 7.  of  the  reading. 

The  5-hole  prism  probe  used  in  this  study  was  calibrated  in  an  open 
jet  wind  tunnel  in  an  effort  to  reduce  errors.  Details  of  the  calibration 
procedure  are  given  in  Reference  ]_  32 / . Treaster  and  Yocum  ]_  58  / con- 

ducted a detailed  statistical  analysis  of  their  data  based  on  a typical 
probability  curve  which  indicated  that  the  50%  error  of  the  mean  was  less 
than  1%  of  the  maximum  value  of  each  probe  coefficient. 

The  inlet  velocity  profiles  were  important  in  the  calculation  pro- 
cedure of  the  secondary  flows  and  therefore  a large  number  of  data  points 
were  obtained.  The  profiles  were  calculated  from  total  pressure  tubes 
and  static  pressure  probes  positioned  in  front  of  the  rotor.  In  all, 
three  total  pressure  rakes  and  one  static  pressure  rake  were  used.  As  the 
rakes  were  rotated,  pressures  were  recorded  at  2°  increments  for  a total 
of  360°.  The  static  pressure  distribution  through  the  boundary  layer  was 
established  by  matching  the  static  pressure  probe  data  to  the  static 
pressure  measurements  on  the  surface.  Each  inlet  velocity  data  point  was 
an  average  of  180  points.  The  accuracy  of  the  measurements  is  determined 
by  the  accuracy  of  the  United  Sensor  static  pressure  probe  which  is  1%  of 
the  reading. 

Desinent  cavitation  numbers  for  the  vortex  were  calculated  from 


pressures  obtained  from  a pitot-static  probe.  The  probe  was  manufactured 
by  United  Sensor  and  the  accuracy  was  1%  of  the  reading.  However,  there 
is  always  data  scatter  for  cavitation  data  because  a desinent  cavitation 
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point  is  visually  determined  and  is  influenced  by  fluid  property 
variations.  The  maximum  scatter  of  the  cavitation  data  is  10%  of  mean 


value. 
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